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ABSTRACT 


Microwave  attenuation  technique*  are  used  to  find  the  combustion 
chamber  temperature  of  a solid  rocket  motor,  and  the  predicted  value 
calculated  by  a theoretical  specific  Impulse  computer  program  and  the 
value  calculated  with  the  microwave  measurements  agree  within  three 
percent.  The  propellants  used  In  this  Investigation  were  seeded  with 
five  percent  potassium  perchlorate  by  weight  to  furnish  enough  electrons 
In  the  plasma  to  cauae  measurable  attenuation.  The  propellants  were 
also  seeded  with  5 to  10  percent  aluminum  to  produce  suspended  metallic 
particles  In  the  plasma. 

The  temperature  profile  inside  the  combustion  chamber  is  deter- 
mined for  a distance  of  two  Inches  away  from  the  burning  surface  by 
scanning  tne  focused  microwave  signal  transverse  to  the  combustion 
chamber  and  by  utilizing  the  motion  of  the  burning  surface  with  respect 
to  the  sntenna3.  The  experimental  and  theoretical  results  for  the 
temperature  gradient  have  a maximum  angular  difference  of  25  degrees 
for  5 percent  aluminized  propellant  and  18.5  degrees  for  10  percent 
aluminized  propellant. 

The  acouscic  oscillation  Inside  a cylindrical  end  burner  is 
theoretically  derived  and  experimentally  observed.  The  oscillation 
frequencies  observed  range  from  3.2  to  4,4  kHz,  whereas  the  theoretical 
oscillation  frequencies  range  from  2.98  to  5.13  kHz  for  various  oscil- 
lation modes. 

Acoustic  gain  and  loss  expressions  are  derived  and  applied  to 
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the  rocket  firings,  The  results  show  that  for  s stable  system,  the 

r 

acoustic  loss  exceeds  the  acoustic  gain  for  some  value  of  A^»  the  real 
component  of  burning  surface  admittance,  and  for  an  unstable  system, 
the  acoustic  gsin  exceeds  the  loss  for  the  same  value  of  A^. 
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I.  INTRODUCTION 


Temperatures  of  a hot,  multispecies,  confined  plasma  are 
difficult  to  measure  with  thermocouples  and  other  pyrometers.  The 
idea  of  determining  the  temperature  of  hot  gases  from  their  radiation 
is  at  least  as  old  as  the  century.  Optical  methods  of  gas  pyrometers, 
using  visible1  and  infrared2'3  radiation,  were  used  in  the  early  1900’s. 
More  complex  spectroscopic  methods  for  gas  temperature  measurement  were 
developed  with  the  advance  in  quantum  theory  and  its  interpretation  of 
atomic  and  molecular  spectra.  Much  literature  is  available  on  these 
methods. *4-7  The  usefulness  of  these  spectroscopic  methods  heavily 
depends  upon  the  understanding  of  what  to  measure  and  how  to  interpret 
the  results.  This  is  mainly  due  to  the  complexity  of  gas  radiation.11 

In  1964,  Jones,  Johnson,  and  Grow9  developed  a microwave  attenua- 
tion measurement  technique  which  was  applied  to  the  plume  of  rocket 
exhausts.  Propagating  microwave  signals  through  the  plasma  and  measur- 
ing the  wave  attenuation  and  phase  shift  provide  information  about  the 
plasma  temperature,  pressure,  and  electron  density.  Olsen  and  Grow10 
further  applied  this  technique  to  the  study  of  temperature  inside  the 
rocket  combustion  charaher  for  the  burning  of  solid  propellant  of  various 
electron  concentrations. 

Metallic  particles  are  often  added  to  solid  propellant  to  in- 
crease the  specific  impulse  of  the  engine  or  to  stabilize  the  combustion 
process  or  both.  The  presence  of  metal  additives  produces  graybody 
radiation  from  the  particle  cloud  in  the  flame  in  addition  tc  the  spectral 
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bands.  These  properties  enable  the  fast  spectral-scanning  pyrometer  to 
measure  the  temperature  of  solid  propellant  flame  containing  suspended 
solid  particles  by  the  two-color  method,11  and  Tour in 1:1  has  noted  that 
only  a aooty  flame  is  sufficiently  gray  for  the  two-color  method  to  be 
applicable.,  Tourin12  alao  noted  that  most  gasea  are  not  gray. 

Part  of  the  present  investigation  is  to  study  the  applicability 
cf  using  m;crowave  measurement  techniques  to  determine  the  combustion 
chamber  temperature  of  solid  propellant  with  various  amounts  of 
aluminum  seeding.  By  measuring  the  chamber  pressure  and  wave  attenuation 
that  is  caused  by  \ known  number  of  free  electrons  in  the  plasma,  the 
combustion  chamber  temperature  of  the  rocket  motor  is  determined  ac- 
curately. 

The  localized  plasma  temperature  can  be  measured  by  focusing  the 
microwave  signal  to  a few  millimeter-beam  radius  and  scanning  it  acroaa 
the  plasma.  By  using  this  method  the  temperature  change  away  from  the 
propellant  burning  aurface  i,s  determined  for  a distance  less  than  two 
inches. 

Acoustic  oscillation  inside  the  rocket  combustion  chamber  is  a 
phenomenon  associated  with  pressure  oscillation  at  the  frequencies 
corresponding  to  the  acouatic  frequencies  of  the  chamber.  In  this  in- 
vestigation, the  acoustic  oscillation  inside  a cylindrical  end  burner 
is  explored  both  theoretically  and  experimentally.  To  the  author’s 
knowledge,  the  oscillation  phenomenon  has  never  been  observed  for  the 
cylindrical  end  burner  rocket  motor. 

In  Chapter  II,  the  theoretical  results  obtained  for  electro- 
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magnetic  wave  impinging  on  a single  electron  source  plasma  will  be 
discussed.  The  dispersion  equation  describing  the  propagation  of  the 
wave  through  the  plasma  region  is  first  derived.  These  mathematical 
relations  -ill  describe  the  wave  attenuation  in  terms  of  the  frequency 
of  the  wave,  the  electron  density  of  the  plasma  and  the  electron  drift 
velocity. 

In  Chapter  111,  the  oscillation  of  acoustic  waves  is  discussed. 
From  the  equation  of  conservation  of  energy,  conservation  of  momentum, 
and  equation  of  continuity,  the  oscillatory  pressure  wave  equation  is 
derived.  The  oscillation  frequency  as  well  as  acoustic  gain  and  loss 
mechanism  are  also  discussed. 

In  Chapter  IV,  the  heat  loss  inside  the  combustion  chamber  is 
discussed  from  the  thermodynamics  aspect.  The  results  are  used  as  a 
comparison  for  the  temperature  profile  measurement. 

All  the  theories  discussed  prev  ously  are  applied  to  the  rocket 
motor  designed  for  this  investigation  in  Chapter  V.  The  experimental- 
setup  and  procedures  are  discussed  in  detail  and  the  theoretical  results 
are  also  presented. 

Chapter  VI  contains  all  the  experimental  Jesuits  and  their 
interpretations . 

The  summary  and  conclusions  are  given  in  Chapter  Vll. 
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II.  ELECTROMAGNETIC  WAVE  PROPAGATION  THROUGH 
A HOT  MULTI SPECIES  PLASMA 


In  this  section,  the  phenomenon  of  microwave  Interaction  with 
a high  temperature,  multlspecles  plasma,  generated  by  the  burning  of 
solid  rocket  propellant.  Is  Investigated  to  determine  the  attenuation 
of  the  Incident  wave.  The  dispersion  equation  Is  derived  from  Max- 
well's equation  and  Langevln's  force  equation,  and  the  attenuation  Is 
found  in  terms  of  the  frequency,  electron  density,  and  drift  velocity. 
The  mean  electron  drift  velocity  is  shown  to  be  a function  of  the 
species  of  the  gas  present,  their  electron  collision  cross  section, 
gas  pressure,  temperature,  and  frequency. 

The  propellant  used  in  this  study  contains  five  per  cent 
potassium  perchlorate  to  furnish  enough  electrons  to  swamp  the  random 
generation  of  free  electrons  and  produce  a measurable  attenuation  in 
the  pressure  range  of  interest.  Metallic  particles  are  very  often 
added  to  propellant  either  to  increase  the  specific  impulse  of  the 
engine  or  to  stabilize  the  combustion  process  or  both.  Various  amounts 
of  aluminum  were  "seeded"; to  produce  a "metallic  cloud"  in  the  combus- 
tion chamber  for  the  purpose  of  this  study.  The  final  form  of  the 
sttenuatlon  expression  is  expressed  in  terms  of  the  wave  frequency,  the 
composition  of  the  gas,  electron  density,  gas  density,  and  temperature. 
Part  of  the  following  derivations  were  given  by  Jones,  Johnson  and  Grow. 

2.1  Dispersion  Equation 

The  propagation  of  electromagnetic  wave  through  a confined 
conducting  medium  is  governed  by  Maxwell's  equation,13 


V 


7 x H 


(2.1) 


JT  + D 


V x E **  -B 


(2.?' 


where  the  total  electronic  current  density,  JT,  is  defined  by 


JT  ■ Pv 


(2.3) 


The  velocity,  v,  in  Eq.  2.3  is  the  average  drift  velocity  of  the  elec- 
tron caused  by  the  impressed  field  and  is  not  the  random  thermal 
velocity  of  the  medium.  From  the  relations  above,  a wave  equation  may 
be  developed  for  the  electric  field  intensity  given  by 

7^E  “ upv  + peE  (2.4) 


This  is  a simple  scalar  equation  which  can  be  used  from  the  geometry 
of  the  test  rocket  motor  and  the  frequency  of  the  applied  signal.  Look- 
ing ahead  momentarily  at  this  point.  Fig.  2.1  shows  the  relative  posi- 
tion of  the  horn  antenna,  dielectric  lens,  and  rocket  combustion  chamber. 

The  operating  frequency  of  the  incoming  signal  is  23.6  GHz,  which 
gives  a free-space  wavelength  of  0.5004  inch.  The  quartz  window  aperture 

measured  one  by  two  inches.  Although  the  plane  wave  approximation  would 
2 

be  valid  for  2L  A ■ 50  inches  for  a 2.5  inch  aperture,  the  directivity 
of  the  horns  is  sufficiently  high  that  transmission  between  the  two  lens- 
fitted  horns  is  only  0.6  dB  less  than  solid  waveguide  and  transmission 
after  insertion  of  the  motor  with  quartz  windows  only  decreases  by  an 
additional  1.0  dB  (see  Appendix  III).  These  results  give  confidence  in 
the  plane  wave  assumption  used  in  the  theory. 
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Relative  position  of  horn  antenna,  dielectric  lens,  and  rocket  combustion  chamber. 


To  find  the  average  drifting  velocity,  v,  in  a sufficiently  high 
density  gaa  molecule  plasma,  the  Langevin  force  equation  is  used  with  a 
damping  term  to  account  for  the  energy  loss  due  to  collision. 14  Since 
the  collision  is  random,  it  is  necessary  to  write  this  equation  for  the 
average  electron;  that  is,  one  thst  experiences  one  collision  per  mean- 
free-path  length.  Thus  we  have 


v + gv  « - E 
* m 


(2.5) 


which,  for  a given  harmonic  time  dependence,  has  a solution  of 


e E 

v *= — 

B g + ill) 

If  the  electric  field  has  a one— dimensional 


(2.6) 

space  dependence  given  by 


E - 


(2.7) 


then  the  simultaneous  solution  of  the  Langevin  equation,  Eq.  2.6,  the 
wave  equation,  and  Eq.  2.7  for  the  propagation  constant  T,  yields  the 
diapers ion  equation, 


r2  + k2  - k2  lw 


‘p  g + iw 


(2.8) 


where  k is  the  propagation  constant  in  free  space  and  kp  is  defined 

by 


u 

k 

P C 


(2.9) 


The  quantity  is  the  plasma  frequency  and  ie  given  by 


(2.10) 


P 


The  quantity  p Is  the  electron  charge  density  of  the  medium  In 

G 

Coulomb  per  unit  volume,  with  the  usual  Identification 

r = a + 16  (2.11) 


where  a is  the  attenuation  constant  in  nepers  per  meter  and  6 is  the 
phase  constant  In  radians  per  meter. 

The  parameter  g in  the  Langevln  equation  has  the  dimension  of 
frequency.  For  a simple  plasma  where  g is  independent  of  electron15 
velocity,  it  reduces  to  the  collision  frequency  of  the  electrons. 

The  rocket  combustion  chamber  plasma  assumed  constitutes  a very  hot, 
dense  mixture  of  gases  and,  under  this  condition,  the  quantity  g may 
be  complex  and  unequal  to  the  collision  frequency.  In  order  to 
allow  for  this  possibility,  the  solution  to  Langevln' s equation  may 
be  written  as 


v - — (B  - ID)  (2.12) 

in 

using  Molmud's  notation.15  The  conductivity,  o,  is  defined  by  the 
equation 


J = oE  (2.13) 

Combining  Eqs.  2.3,  2.10,  2.12,  snd  2.13  yields  Eq.  2.14, 
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(2.14) 


o ■ ui*  e (B  - ID) 

P ° 

From  Eq.  2.14  it  is  clearly  seen  that  "B"  and  "D"  are  directly  related 
to  the  conductivity  of  the  plasma.  Using  this  relation  the  resultant 
determinantal  equation  is  given  by 


T2  + k2  - ik2wB  - k^u)D  - 0 
P P 


(2.15) 


Substituting  F from  Eq.  2.11  and  collecting  real  and  imaginary  parts 
c‘.  Eq.  2.15,  we  have 


f 


a 


<i) 


(2.16) 


6 


U) 


(2.17) 


From  these  equations  a and  B are  functions  of  wave  frequency, 
plasma  frequency,  which  is  related  to  electron  density,  and  the  complex 
drift  velocity  terms,  B and  D. 


2.2  Investigation  of  the  Drift  Velocity 

In  order  to  complete  the  derivation  of  attenuation  and  phase 
constants  in  terms  of  physical  parameters,  the  functional  dependence 
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of  drift  velocity  must  be  investigated.  The  velocity  term  described 
in  Gq.  2.12  is  an  average  velocity  of  all  electrons  in  the  gas  and  is 
in  the  same  direction  aa  the  electric  field.  If  a coordinate  system 


is  assumed  so  that  the  electric  field  is  in  the  z-direction,  then  the 
velocity  given  by  Eq.  2.12  is  actually  the  average  velocity  in  the  z- 
direction.  Using  Maxwell's  law  of  the  velocity  distribution,  the 
average  drift  velocity  is  written  as 


| V2fd, 

o 

CD 

| fdT 
O 


(2.18) 


where  the  function  f must  satisfy  the  Boltzmann  equation  which  is  a 
mathematical  relation  that  specifies  the  phase-space  distributl  i 
function  for  a particular  system  of  particles.17 


3f 

3t 


(2.19) 


MargenaulH  has  shown  that  the  solution  to  Eq.  2.19  has,  as  the  first 
terms  of  its  expansion, 


f = fQ  + ~ v2  cos  wt  + sin  wtj 


(2.20) 


where  fQ,  f^,  and  g^  are  coefficients  of  the  expansion.  Equation  2.18 
may  then  be  rewritten  as 
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2ir  n z 


| v^  ^f^  cos  wt  + g^  sin  <Dt^  vz  v 

t>»0  0-0  vO J_ 

2n  n “> 

| | fQ  v^dv  sin  0d0d$ 


‘dv  sin  Od'td'l 


4>=0  e-o  v-0 


(2.21) 


Relations  may  be  developed16  between  f f^,  and  g^  such  that  Eq.  2.21 
may  be  expressed  in  complex  form  as 


2 ii  « 


f f [ eE  4,  Amv  f m \ ' (,  iAw\  2„  . 

T V dV  2 2,2  n\2krJ  l1  - — j COS  6 6l"  0d’d* 

~ J-  ~ V + uj  A v ' ' ' 


■0  0=0  v=0 


2"  " z 


f f "fe) 


2kT  2 

e v sin  0d0d$dv 


<p“0  0*0  v=0 


(2.22) 


where  A is  the  mean  free  path  between  collisions,  and  n is  the  electron 
density.  If  a change  of  variable  is  made  such  that 


v 

y "7 


(2.23) 


and  the  0 and  $ integrations  are  carried  out,  then  Eq.  2.22  becomes 


8n  eE  ( / B \3/2  vA  e~Bv 

3 n ) l«/  iui  + v 


(2.24) 


where 


(2.25) 
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For  most  cases  v is  a velocity  dependent  function.  For  ao-called 

Maxwellian  gas,  in  which  the  potential  between  electrona  and  neutral 

-4 

molecules  varied  as  r , v is  not  a velocity  dependent  function. 

A mixture  of  frequency  is  found  from  the  average  electron-molecule 
collision  cross  section,  Qave«  given  by 

Q - }'  m.  Q.  (2.26) 

ave  L i M v 

where  the  m^'s  are  the  mole  fraction  averaging  factors  of  each  species 
of  gas  present,  since  the  collision  frequency  of  the  electrons  and 
a particular  neutral  gas  is  given  by 

v1  = ovQ1  (2.27) 


where  p is  the  gas  particle  density.  The  average  collision  frequency 
for  a given  electron  velocity  is 


v 

avg 


(2.28) 


Substitution  of  Eq.  26  in  Eq.  28  gives 


v 

avg 


pv  l mi 


(2.29) 


The  average  drift  velocity  of  electrons  due  to  the  incident 
microwave  signal  as  given  in  Eq.  2.24  is  now  a function  of  the  fre- 
quency of  the  wave,  the  mole  fraction  of  each  gas  species,  and  each 
collision  cross  section.  Rewriting  Eq.  2.24  we  get 
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4 -6v  . 

v e dv 

iu>  + pv  £ Qj 


(2.30) 


To  obtain  useable  expressions  for  the  cross-section  terns  In 
the  summation  In  Eq.  2.30,  a power  series  approximation  of  the  col- 
lision cross  sections  will  be  used.  This  approximation  Is  valid  for 
only  a limited  range  of  velocity;  however,  It  Is  known  that  the  varia- 
tion In  the  combustion  gas  temperature  for  pressures  between  20  psla 
and  100  psla  is  only  about  1008  K out  of  an  equilibrium  temperature  of 
2600°  K,  Indicating  that  the  thermal  energy  of  the  electrons  actually 
changes  very  little  over  the  entire  pressure  range  considered.  It 
will  be  shown  In  a later  section  that  the  significant  collision 
species  In  the  plasma  are  water,  hydrogen  chloride,  carbon  dioxide, 
carbon  monoxide,  diatomic  hydrogen,  and  diatomic  nitrogen.  Abundant 
Information  exists  about  the  dependence  of  the  collision  cvoss  sec- 
don  on  electron  energy  for  these  substances.  Using  the  fact  that 

the  thermal  energy  of  the  plajma  electrons  is  of  the  order  of  0.22 
electron  volts  and  varies  only  a few  percent  from  that  value,  the 
following  expressions  represent  good  approximations  to  the  collision 
cross  sections  for  the  given  gaseous  species: 


5.90  x IQ 

2 

v 


(2.31) 


n » 1.835  x 10 
"HC1  2 


(2.32) 
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(2.33) 


Qrn  - 7.16  x 10 
C02  v 


-14 


-20  ( —25  ^ 

Qco  » 2.65  x 10  + {,2.08  x 10  jv  (2.34) 


Qp  * 8.92  x io-20  + (l.464  x lo“25)  v (2.35) 


Q = (2.70  x io'25)  v 
2 


(2.36) 


where  v is  the  thermal  velocity  of  the  plasma  electrons.  It  is  seen 
that  the  summation  in  Eq.  2.29  contains  terms  of  the  form 


Qv  = const.  =•  c^ 


(2.37) 


where  the  value  of  r is  2,  1,  0,  and  -1.  Equation  2.30,  when  written 
so  that  the  constants  of  Eq.  2.37  appear,  becomes 

® 2 

- . 8^eE  . fly2  ( v5  e"gV  dv 

3 ni  J . . 2.  3 

io)v  + m0c0p  + m,  c,ov  + m c pv  + m,c,pv 
O Z Z 11  OO  —1  —1 


(2.38) 


Making  the  change  of  variable  such  that 


y - 6 1/1 2 v 


(2.39) 


and  substituting  Eqs.  2.39  and  2.12  into  Eq.  2.38  gives 
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co 


B - ID 


8 

3A  j 


5 -y 

By  e y dy 


3/2  1/2  2 3 

iuiBy  + m_c0pB  + ra.c.pBy  + m c pB  y + m .c  .py 

L L 11  00  -1-1 


(2.40) 


Upon  rationalizing  the  denominator,  it  is  seen  that 


B = 


“ { a 3/2  . „ . 1/2  2 3\  „ 5 -y2  . 

8 [ (m2c2pS  + miclp6y  + moCopS  y + m^c^py  j By  e dy 

J 2 l 3/2  1/2  2 3\2 

*>  (u)By)  + ^m2c2pB  + raicjpBy  + “0C0PB  y + m_1c_1py  ) 

(2.41) 


and 


3/fi 


fl2  6 -y  . 

tiB  y e J dy 


o (uj  By ) 


("2c2f’ 


3/2 

B + ra.c.pBy  + m c PB 

11  o o 


1/2  2 . 3 \ 

y + m_1c_i°y  / 


(2.42) 


With  the  determination  of  B and  D as  functions  of  the  incident 
wave  frequency,  gas  composition,  and  temperature,  the  propagation  con- 
stants of  Eqs.  2.16  and  2.17  are  determined  except  for  thu  ”alue  of 
the  plasma  frequency.  Once  the  electron  density  is  known  for  the 
plasma  of  a given  propellant,  the  wave  attenuation  is  completely 
determined  and  can  be  theoretically  predicted  for  a given  geometry. 

2.3  Investigation  of  the  Electron  Density 

In  this  derivation  of  the  propagation  constants  for  the  incident 
electromagnetic  wave,  the  temperature  contained  in  the  expressions 
thus  far  is  the  temperature  of  the  electrons  and  was  introduced  through 


‘-'k,-  .Si,  • 
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the  mean  velocity  expression,  Eq.  2.12.  When  the  electrons  are  in 
thermal  equilibrium  with  the  gas  molecules,  the  temperature  in  a and 
3 can  be  equated  to  the  gas  temperature,  which  in  turn  relates  to 
the  gas  density  through  the  perfect  gas  law, 


P * 


(2.43) 


where  Nq  is  Avogadro's  number,  R is  the  gas  constant,  and  p is  the 
pressure.  The  perfect  gas  law  then  allows  the  expression  of  density 
in  terms  of  temperature  provided  the  pressure  is  known.  The  equilibrium 
electron  density  then  may  be  derived  from  the  thermal  state  of  the  gas. 
It  will  be  assumed  that  the  source  of  electrons  is  the  seeded  potas- 
sium perchlorate  and  that  the  mechanism  of  free  electron  production 
is  thermal  ionization.  There  are  two  ways  of  considering  the  produc- 
tion of  electrons;  they  are: 

1.  Assuming  that  no  sinks  for  electrons  exist  in  the  gas. 

2.  Assuming  that  chlorine  acts  as  a sink  for  electrons. 

2.3.1  Assuming  No  Electron  Sinks 
For  the  simple  reaction 

K + K*  + e"  (2.44a) 

A1  ->■  Al+ + e"  (2^4b) 

Saha's  equatior^1  may  be  written  as 
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where  N is  the  total  number  density  per  unit  volume  of  potassium  and 
aluminum,  X is  the  fraction  of  those  atoms  ionized,  k is  Boltzmann’s 
constant,  h is  Plank's  constant,  mg  is  the  electronic  mass,  e is  the 
electronic  charge,  and  I is  the  ionization  potential.  For  a given 
temperature,  the  number  density  for  electrons  per  unit  volume  is 
simply 

ne  = XN  (2.46) 


and  the  plasma  frequency  defined  in  Eq.  2.10  becomes 


2 

(0 

P 


(2.47) 


2.3.2  Consideration  of  Electron  Sinks 

When  electron  recombination  is  taken  into  account,  not  only  the 
number  density  of  the  electron  source,  but  also  the  number  density  of 
the  electron  sink  must  be  known.  Chlorine  will  be  assumed  to  be  the 
only  effective  electron  sink  because  it  is  the  only  constituent  present 
in  substantial  amounts  in  the  propellant  that  has  an  electron  affinity 
of  the  order  of  the  low  ionization  potential  of  potassium. 

With  the  following  two  reactions  considered,  the  chlorine  elec- 
tron sink  is  defined  as 
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■ ■*  .wttMjltt 


k 


m 


K -*■  K+  + e 


(2.48a) 


A1  ■>  Al+  + e"  (2.48b) 

Cl  + e~  -»■  Cl"  (2.48c) 

The  chemical  equilibrium  constants  for  this  pair  are  given  by 


.Al+][e~]  [k+]  [e~] 

(Al]  + [K] 


(2.49a) 


[Cl3  iel  =a 
[a  ] ' 2 


(2.49b) 


When  the  bracketed  quantities  denote  concentrations  in  number  per  unit 
volume,  the  equilibrium  constants  may  be  written  as 


(2.50a) 


and 


a 


2 


= 4 


(2.50b) 


where  1^  and  I are  the  ionization  potential  of  potassium  and 
aluminum,  respectively,  and  is  the  electron  affinity  of  chlorine.2^ 
Defining  the  total  quantities  of  the  reacting  elements  by  the 
following  equations: 
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A 


1 


(2.51a) 


and 


A2  = [Cl]  + [ci“]  (2.51b) 

it  is  possible  to  combine  Eqs.  2.49  and  2.51  with  the  neutrality  con- 

dition, 

[e  ] + [Cl]  = [k+]  + [All  (2.52) 

to  form  a cubic  equation  in  electron  density, 

n3  + An2  + Bn  + C * 0 (2.53) 

where 

n = [e~J  (2.54) 

and  the  coeff ici«,«*’s  are  given  by 

A = + a2  + A2  (2.55a) 

B « a^A2  - A:  + o2)  (2.55b) 

C = -°'1a2Ai  (2.55c) 


The  determined  value  then  establishes  the  plasma  frequency  through 
Eq.  2.40. 
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With  the  presence  of  aluminum  in  the  solid  propellant,  consider- 
able amounts  of  aluminum  oxide  particles  are  present  in  the  plasma. 

The  amount  of  A^O^  at  various  conditions  is  listed  in  Appendix  I. 

Each  aluminum  oxide  particle  in  the  plasma  can  be  considered 
as  an  isolated  positive-charged  particle  since  the  average  separation 
distance  is  about  7 microns,  as  reported  by  Richardson.24  Compared 
to  the  1 micron  particle  radius,25  the  interaction  particles  can  be 
neglected.  Each  isolated  particle  is  surrounded  by  a thin  sheath  of 
electrons  with  a potential  <}>  as  a function  of  distance  r given  by 
Jones. 1 


t = B 


2p  e 
o_ 

c kT 
o 


1/2 


(2.56) 


where  B is  an  arbitrary  constant  and  Pq  is  the  charge  density. 

An  exact  determination  of  the  charge  density  profile  has  been 
reported  by  Johnson  and  Bullock.  They  concluded  that  the 
electrons  contained  in  the  thin  sheath  are  not  contributing  to  the 
average  electron  density  in  the  plasma.  Thus,  this  effect  can  be 
neglected  in  the  calculation  of  microwave  attenuation  when  propagating 
through  the  plasma. 

2. A Plasma  Temperature  from  the  Attenuation  Constant 

To  relate  the  attenuation  of  a microwave  beam  In  a plasma  to 
the  temperature,  consider  that  a plane  wave  normally  impinges  upon  a 
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plasma  region  and  passes  through  the  plasma  Into  free  space  as  Il- 
lustrated in  Fig.  2.2.  At  plane  I In  the  figure,  due  to  differences 
in  the  characteristic  impedance  of  the  plasma  and  free  space,  a re- 
flection  coefficient  may  be  defined  in  terms  of  the  two  impedance-  by 

> 

* ■ rrnr  <2-57> 

1 o 

where  Zq  is  the  characteristic  impedance  of  free  space  and  Zj  is  that 
of  the  plasma.  If  it  is  assumed  that  the  reflection  coefficient  is 
small  compared  to  unity,  then  the  amplitude  of  the  transmitted  wave, 
an  infinitesimal  distance  to  the  right  of  plane  I,  will  be  essentially 
equal  to  the  amplitude  of  the  free  space  wave.  If  the  transmitted 
amplitude  of  the  electric  field  is  designated  as  Ep  and  the  incident 
wave  by  Eq,  then  the  amplitude  of  the  electric  field  within  the  plasma 
is  given  by 

E - E e"Fz  (2.58) 

p o 

where  z is  measured  to  the  right  of  the  plasma  boundary  I.  If  the 
assumption  of  a small  reflection  coefficient  ia  again  made  at  plane 
II,  then  the  electric  field  amplitude  at  plane  II  designated  by  E^ 
is  given  by 

Et  “ Eq  e_rL  (2.59) 

The  attenuation  of  the  transmitted  wave  in  decibels  by  definition  is 
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Fig.  2.2.  Diagram  of  the  assumed  electromagnetic  field  geometry. 


(2.60a) 


A 


20  108io  r 

o 


or 


A 


20  log1Q  e 


-PL 


(2.60b) 


In  terma  of  the  magnitude  of  a of  Eq.  2.16, 

A ■ 20  aL  log^Q  e 


(2.61) 


Substitution  of  Eq.  2.16  Into  Eq.  2.61  yielda 


A = 10/2  kL  log1Q  e 


(2.62) 


It  la  seen  that  one  attenuation  meaaurement  through  a specified  path 
length  of  plaama  of  a known  compoaition  and  pressure  can  yield  the 
plasma  temperature. 
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III.  PRESSURE  WAVE  OSCILLATION  INSIDE  THE  COMBUSTION  CHAMBER 


The  equation  describing  pressure  wave  or  acoustic  wave  oscilla- 
tion Inside  the  combustion  chamber  of  a rocket  motor  la  derived  from 
the  continuity  equation  and  the  energy  and  momentum  conservation  equa- 
tions. Several  valid  and  practical  aaaumptlona  are  made  to  simplify 
the  analysis.  The  equation  for  the  oscillation  it  then  solved  by 
Creen's  function  and  perturbation  methods..  The  final  equations  for 
both  traveling  and  standing  waves  are  obtained  In  the  lowest  order  of 
gas  mean  flow  Mach  number  of  cylindrical  geometry. 

The  growth  rate  of  acoustic  waves  is  expreaaed  in  terms  of  the 
pressure  wave,  the  flow  velocity  of  the  gas,  and  the  acoustic  admit- 
tance of  the  burning  surface.  The  losa  mechanism  involved  not  only  the 
nozzle  size  but  also  the  damping  effecta  of  solid  particles. 

Part  of  the  derivationa  can  be  found  in  Culick'a  work. 27 


3. 1 Ceneral  Equation 

The  motion  of  a viscous,  compressible,  heat-conducting  gas,  as 
in  the  combustion  chamber  of  a rocket  motor,  is  given  by  the  following 
equations: 

Continuity  equation: 

- 77  + V • (on)  - 0 (3.1) 


Conservation  of  momentum  or  Navier-Stokes  equation: 


P 


du  VP 
dt  y 


F + F. 
s b 


(3.2) 
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Conservation  of  energy: 


p dF“ 


_ ^ ^jp 

— — ^ * a function  of  quantities  that  relate  fluid 
viscosity  to  quantities  such  as  fluid  flow 
velocity,  specific  heat,  thermal  conductivity, 


and  the  boundary  layer  parameter. 


(3.3) 


Equation  of  state: 


P = PT 


(3.4) 


where  the  notations  used  are  listed  as  follows: 

F * viscous  force 
s 

= body  force,  related  to  the  particle  damping  effects  on 
the  growth  of  the  acoustic  wave 
p * density  of  the  gas 
u * velocity  vector 
Y ■ ratio  of  specific  heat 
P “ amplitude  of  pressure 
T = temperature 

Equations  3.1  through  3.4  are  written  in  terms  of  the  dimensionless 
variables;  i.e., 


_ P’ 

P . r 
o 


u = — 
— a 


(3.5) 
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aQ  is  the  average  speed  of  sound,  all  "prime”  superscripts  refer  to 
dimensional  quantities,  all  "zero"  subscripts  refer  to  average  (dimen- 
sional) conditions  in  the  absence  of  waves. 

With  several  crucial  and  valid  assumptions,  the  work  will  be 
greatly  simplified: 

1.  The  linear  burning  rate  is  always  small  compared  to  the 
wave  propagation  speed.  Generally  speaking,  the  propel- 
lant useo  in  the  experiments  always  burned  less  than  half 
an  inch  per  second,  whereas  the  speed  cf  sound  is  around 
3000  feet  per  second.  Hence,  the  rate  of  recession  of  the 
surface  can  be  ignored  in  a treatment  of  steady  waves. 

2.  The  combustion  occurs  within  a very  thin  region  on  the 

-1  -3 

order  of  10  to  10  inches  thick.  So  the  coupling 
between  the  oscillations  and  the  burning  may  be  safely 
introduced  as  a boundary  condition  placed  on  the  waves. 

3.  Mach  number  representative  of  the  mean  gas  flow  within 

the  chamber  at  the  burning  surface  is  small;2 / i.e.,  the 

Mach  number  at  burning  surface  is  • 0.05  and  the  Mach 

number  at  exit  is  M„  * 0.2  to  0.3. 

e 

Before  using  the  perturbation  method  to  solve  these  equations, 
we  must  also  assume  that  the  average  Mach  number  of  the  acoustic 
disturbance,  c,  is  small,  and  the  Mach  number  of  the  mean  flow  field 
at  the  burning  surface  where 

V - (r)  <3-6> 

' chamber  wall 
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is  also  small. 


The  mean  flow  within  the  chamber  ia  represented  by 

(u)  mean  flow  - vU  (3.7) 

where  U reflects  the  geometrical  aspect  of  the  mean  flow  pattern. 

The  effect  of  viscosity  could  be  aasumed  to  be  confined  to  a 
thin  layer  of  chamber  walls,  and  the  viacoua  terms  are  therefore 
dropped  out  from  the  governing  equation.  Thus  we  get 

+ V • (pu)  = 0 (3.8) 


P 


3“  VP 

at  v 


(3.9) 


P 


— - rjil  dp  „ 

dt  V dt  u 


(3.10) 


Differentiating  Eq.  3.4  and  combining  with  Eq.  3.10,  we  obtain 


1 dP  _ dp 
y P " p 


(3.11) 


v 

Hence  we  have  P ■ o which  indlcatea  that  the  inviscid  part  of  the 
flow  may  be  regarded  as  iseniropic.  Using  this  result,  Eq.  3.8  can 
be  written  as 


— + U • VP  + yPV 


(3.12) 
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(3.13) 


These  equations  ate  the  governing  equations  ,or  the  inviscld 
chamber  motion. 

To  apply  the  perturbation  nethod,  the  folloving  perturbation 
expansions  are  used: 


P-l  + tV(1>  + E2p«>  + ... 
T ■ 1 + tT(l)  + c2T«)  t ... 


(3.14  a) 
(3.14b) 
(3.14c) 


u >=  vU  + eu(1)  + e2u(2)  + . . 

Substituting  these  values  into  Eqs.  3.12  and  3.9,  „e  get 
Order  of  r: 


(3.14d) 


St 


+ 


(3.15) 


Su^  yp  (1)  r 

— + v[u  . 7u(1)  + u(1)  . vul  + F(1) 

— ■*  b 


(3.16) 
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Order  of  e : 


3P 


+ yV  • u(2)  - -yP(1)7u(1)  - u(1)  • VP*1*  - v[u  • VP*2)] 

(3.17) 


3u(2)  . VP(2)  (1)  ^yL)  (1)  „ (1)  T (2)  „„ 

^r- — + = -P  “Z u • Vu  - V u • vu 

3t  y 3t  — — L—  — 


.(1) 


+ u • Vu*2*  + P^u*1*  • VU  + U • Vu*1*)  + F*2* 


(3.18) 


where 

U is  assumed  to  be  steady  and  solenoidal;  i.e.,  V • U * 0 
Differentiating  Eq.  3.15  with  respect  to  time,  we  obtain 

+ y ~ (v  • u*1*)  - -vfl  • — (vP(1))  (3.19) 

Taking  the  diverge^e  of  Eq.  3.16,  we  get 

yV  • + 7 ’ (7P(1))  - -yv[v  (u  • Vu*1*  + u*1*  • vu)]  + y V • F*1* 

(3.20) 


Combining  these  two  equations,  we  get 


2 (1) 

- vV1*  - -yV 
3t 


F*1*  + yy  V 

b 


- U ■ 


(3.21) 
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In  anticipation  of  oscillation,  we  assume: 


p<i>  . TP<i> 

„<»  - 

F.a)  - «<» 

D — 

where  k is  the  complex  frequency  and 

k = ft  + jA 

The  amplitudes  of  p^\  and  may  be  complex. 

With  these  relations,  Eq.  3.21  can  be  written  as 

V2p(1)  + kV1*  - vga) 

where 

K(1)  - |jk(u  ■ Pp(1>)  + a<»  . ,„)  + (?.) 


Now,  introducing  an  admittance  function,  A, 

A = Af  + jAt 


such  that 


burning  surface 


(3.22a) 

(3.22b) 

(3.22c) 

(3.23) 

(3.24) 


(3.25) 


(3.26) 


(3.27) 
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where  11  a.  > P are  the  complex  velocity  and  pressure  fluctuation  at 
the  wall,  n Is  an  outward  pointing  unit  vector  normal  to  the  chamber 
wall. 

Substituting  Eq.  3.16  we  get 

n • Vp^  = n ♦ yj-dkg/1^  - \( U * * Vlj)  + Xf 

(3.28) 


Also,  using  Eq.  3.27, 


n . q(1>  = -vAp(1) 


vAp  at  boundary 


(3.29) 


where 


at  Inhibited  surfaces 


A^  at  burning  surface28 


— A at  nozzle  exits 

\i  n 


\j  * Mach  number  of  flow  In  nozzle  entrance 
n 


) *' 


A = nozzle  admittance 
n 


So,  the  boundary  condition  can  be  expressed  as 


n . vp(1)  = -vh(1) 


(3.30) 


where 


h(1)  = |-jkAp(1)  + n • y[u  • 7cl(1)  + • VU  - (-)  f (1)]J  (3.31) 
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and  A is  the  appropriate  admittance  function  at  the  chamber  boundary. 

3.  2 Solution  by  Green’s  Function  Method 

Equations  3.24  and  3.30  are  nothing  more  than  nonhomogeneous 
wave  equations  which  satisfy  the  nonhomogeneous  Newmann  boundary  con- 
dition. 

Green’s  function  method  is  used  to  solve  this  Helmholtz  equation; 
the  first  order  solution  of  Green's  function  aatiafles 

(v2  + k2)  G (r|?Q)  = 5^  - ro)  (3.32) 


(->  \ 

r - r j is  the  delta  function  representing  an  harmonic  point 
source  at  the  position  rQ;  thus  G represents  the  disturbance  received 
at  r.  The  boundary  condition  chosen  for  G is  that  the  normal  compo- 
nent of  its  gradient  vanishes  on  the  entire  bounding  surface: 

n • VG  (lIIq)  e 0 for  r * rg  (3.33) 


where  r is  a vector  defining  position  of  polnta  on  the  bounding 
a 

surface. 

Multiplying  Eq.  3.24  by  G ^rjr^j  and  Eq.  3.32  by  p^,  and  sub- 
tracting the  two  equations,  we  get 
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Integrating  over  the  volume  of  the  chamber  — 0<r<l,0<z<  L/R, 
0 < $ < 2n  in  cylindrical  coordinate  where 
R is  the  chamber  radius 

L is  the  characteristic  length.  For  axial  modes  L is  the 
chamber  length;  for  tangential  or  radial  modes,  L is 
the  chamber  radius  or  the  width. 

hence 


v v v 


(3.35) 


The  last  term  is  just  evaluated  at  r = rQ  using  symmetry  property; 

i.e. , 


‘ P'O  ■ [‘  M* 


(3.36) 


where  * denotes  complex  conjugate. 

(1)  -*• 

We  get  p (r)  in  terms  of  G*: 

p(1)(r)  = v | G*  (r|rQ)  g(1)  dv  - | [g*V2p(1;  - p(1)  V2G*]  dv 
v v 

(3.37) 

Using  Green's  theorem  to  convert  volume  integral  to  surface  integral, 
we  have 

j (g*72p(1)  - p(1)V2G*)  dv  - | (c*7p(1)  - p(1)VG*)  • nds 
v s 

(3.38) 
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but 


n 


n • VG*  ■ 0 

over  the  boundary  surface.  Thus 

P(1)(r)  - v j G*  (rjrj  g(1)  dv  + v j G*  (r|rQa)  h(1)  ds  (3.39) 
v a 

► 

where  rQs  refers  to  the  value  at  points  on  the  boundary  surface. 

A useful  expression  for  G^rjr^  ia  the  expansion  in  normal 
modes  of  the  unperturbed  wave  equation.  Assuming 

C - l Aa*a(0  (3.40) 

a 


where 


(’2  + k„2)  ♦„&  ■ 0 


with 

A 

n • V -0  for  r * r 

a s 


and  t stands  for  three  indices:  l,  m,  and  n. 

The  eigenfunctions  must  also  be  orthogonal;  i.e., 

(VM  " f K *6 dv  m o for  a * 6 

v 


(3.41) 


(3.42) 


(3.43) 
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The  normalization  factor  E is 


,2  * 

•-  “ 0 Hi 

a a a 


(3. 44) 


Then  Eq.  3.32  can  be  written  as 


l Aq  (v\  + k2^)  = fi(?  - rj 


(3.45) 


From  Eq.  3.41  we  have 


V2^  (r) 

a 


-k>  (r) 

a a 


(3.46) 


l AJa(h  (k  - k2)  = S[r  - ro) 


(3.47) 


K ~+ 

Multiplying  2 i^ft(r)  by  Eq.  3.47  and  integrating  over  volume  v*  we  have 


l *J<*>  6(r  - <*v  = l *J(r)  1 1 Au*a(r)  (l 


k2  - k2)  dv 


(3.48) 


Using  orthogonal  property  and  by  the  definition  of  6 function,  we  get 

KCO’K  (k2-k«)  fv>  <3-45> 


**  (r  ) = A (k2  - k2)  E2 

n \ oJ  u \ u / a 


(3.50) 
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Therefore 


Hence, 


, . ♦»(*.) 


Wo)'  l 


ffio) 


i E2  U2  - V2) 

a \ a J 


2 * 

l ■ ifi  tj« 
a J a a 


3.3  Solution  to  Zeroth  Order  in  v 


Rewriting  the  governing  equstlons, 


72p<»  < k2p(1)  - vg(1> 


n * = -vh^  at  boundary 


From  these  two  governing  equations,  p csn  be  represented  as  a 
perturbation  expansion  with  v, 


p<1J  - p<10>  + vp(11>  + v2p(12)  + . 


where  the  first  element  of  the  superscript  refers  to  the  order  of  e 
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and  the  second  refers  to  the  order  of  v.  By  the  same  token,  since 
and  h^  are  functions  of  p^,  we  can  also  write 


(3.54) 


Note  Eq.  3.56  is  identical  to  the  boundary  value  problem  shown  by 
Eq.  3-4l.  This  simply  demonstrates  that  the  solution  as  it  i9  un- 
folding is  a perturbation  on  the  classical  acoustic  solution.  Thus, 
. , and  it  is  proper  to  write 

P(1)  = <|/  + 0(v)  (3.60) 

where  0(j)  is  for  higher  order  terms  in  v.  Solution  of  p^  can  be 
found  in  terms  of  Green's  function  .here 
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p(1)  - v | Gg(1)dv  + V j G*h(1)ds 


G-  l 


S E2  (k2  - k2) 
a v a/ 


So, 


(1) 


I (1)  r "a''1  *J;o)  , , f , (1)  , V'>  **('<.) 
"MB  l 2 V 2~~2\  dv  + v h i T71 2\ 

J a E2  (kz  - kz)  i a e:  (V  - K ) 


(*2-V 


*r.(r) 

v L 


‘ Ea  (fc2  - k«) 


(1)  *+  [ 
gU;  *a(r)  dv  + J 


, . (1)  * -k 

(r)  dv  + | h''  t|*  (r)  ds 


For  a particular  mode  of  oscillation,  say  N, 


,(D  _ . , 

N " V 


p 

‘ 

V 

( m * ( (l)  * 

j gU;  *N  dv  + j hU>  ^ ds 

. v s 

2 /,2  , 2\ 

en  Vk  " kJ 

^(r) 

r < i 

+ v l 


Ea  (fc2  - ka) 

To  achieve  Eq.  3.60,  we  need  to  have 


(1)  * (1)  * 

' ; i|>  dv  + hK  ' ds 
o a 


s 


(3.61) 


(3.62) 


ds 


(3.63) 


(3.64) 
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*N  dv  + 


(l)  * 

h(U  *N  ds 


1 


(3.65) 


So , 


k 


2 


j g(1)  dv  + | h(1)  **  ds 

V S 


(3.66) 


As  Eq.  J . 23  shows. 


k = n + JA 


(3.23) 


At  this  point,  we  car  see  Q represents  the  frequency  of  the  fluctua- 
tion, and  A represents  the  rate  of  growth.  Now  if  we  write 


<l0)  + + v2S)(12*  + ... 

(3.67) 

(1°>  + vA(11)  + v2A<12>  + ... 

(3.68) 

expanding  Eq.  3.66  in  Taylor's  series,  we  get 


Si 


2kN^ 


<1)  ...  + f h(D  * 


g-  - 4v  j 


*N  dS 


+ 0 


(v2) 


(3.69) 


and  it  is  apparent  that  k has  no 

Au°) 

Uu°) 


imaginary  part  of  order  unity.  Thus, 


« 0 

(3.70) 

= Si 

(3.71) 
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Also, 


2kNEN 


-<10>  ** 


(3.72) 


A(U)  = — I 6U0)  ^ dv  + f h<10)  C ds  (3.73) 

Mi  .1  l J 


Equations  3.70  and  3.71  show  that  the  frequency  of  the  oscillations 
Is  the  same  as  the  acoustic  frequency  to  zeroth  order  in  the  mean 
flow  Mach  number  v. 

Now,  back  to  Eqs.  3.56  and  3.57;  the  zeroth  order  in  v motion 
is  given  by 


vV10>  + 0<10>  p<10>  - 0 


(10) 

n • pv  "0  at  the  boundary 


In  cylindrical  coordinates,  these  equations  become 


I _L  L 3g(10)>)  + Jl&L 


(10)  2 (10)  2 
2 — + ^ + 0(10)  p(10)  - 0 (3.74) 


on  r * R 


at  z = 0,  z ■ - 


(3.75) 


- 40  - 


Here  it  is  convenient  to  nondimensionalize  the  lengths,  choosing  R as 
characteristic  length,  nondimens ions 1 chamber  radius  and  length  arev 
respectively,  1 snd  L/R.  Then  Eq.  3.74  can  easily  be  solved.  Hence 
w o get 

■ 

J fk  r)  ros  (mO)  cos  (k.z)  standing  wave  (3.76) 

mv  mn  / \ *■  / 

p<10>  - . 

a 

Jm(kmn  r)  6 ^ C0S  (k'z)  traveling  wave  (3.77) 

where 

a * (£,  m,  n) 

+ is  the  sign  for  clockwise  traveling  waves 
- is  the  sign  for  counterclockwise  traveling  waves 
To  satisfy  the  boundary  condition, 

ki - (!)  * = °* 2 ••• 

d J (k  r) 

k is  the  root  ot * - 0 

mn  dr  Jr=1 

For  either  traveling  or  standing  waves,  some  values  are  listed 
in  the  following  table  ss  given  by  Culick.27 

Table  3.1.  Frequency  eigenvalues  k for  cylindrical  chamber. 

mn 

d J ( k r j 

Values  of  the  roots  ■ ™— -J5S = q 

dr 
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The  frequencies  are: 


and  conversion  to  a dimensional  form  gives 


(3.78) 


f' 


mn 


(Hz) 


(3.79) 


The  normalization  constants  for  either  standing  or  traveling  waves 
are: 


m - 0 


m ^ 0 


(3.80) 


With  Eqs.  3.22a,  3.53,  and  3.76,  we  can  easily  write  out  a 
pressure  wave  equation  in  the  form  shown  in  Eq.  3.14d  with  the  growth 
rate  of  acoustic  wave.  A,  yet  to  be  determined. 

3.4  Growth  of  Acoustic  Wave 

In  the  previous  section,  the  pressure  wave  equation  in  cylindrical 
chambers  was  derived.  Rewriting  these  equations,  we  have 

k * a + jA  ’ (3.23) 

and  W represents  the  oscillation  frequency,  A represents  the  growth  nr 
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decay  of  the  oscillatory  wave.  We  also  derived 


51  - 5l<10>  + vSl(11)  + v2fi^12)  + ... 


(3.67) 


A - A(10)  + VA(11)  + vVl2>  + 


(3.68) 


of  which 


A(10)  * 0 


(3.70) 


n(10)  - k„ 


(3.71) 


51(11)  - -U  R 

Vn 


(10) 


(10)  * 


g ' ' dv  + I h ds 


Lv 


(3.72) 


A(11)  - -U  I 

ZkA 


jg(10)*>  + 

1 -<10)  < *• 

_ V 

3 

(3.73) 


with  representing  the  pressure  wave.  Substituting  yp^10^  for 

Pi  N 

and  for  first  order  v,  we  can  write 


k - Sl<10>  + + JA‘U>]  + 0 [\J2] 


(3.81) 


where  o[v2]  represents  higher  order  v terms.  Using  Eqs.  3.72  and  3.73, 


we  obtain 


&(11)  + A(U) 


(R  +J1) 

2Q<10>E2 

a 


f 8(10)  P(10)*dv  + f h(l0>  p(10^ds 


(3.82) 
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and  the  oscillation  fre- 


So  the  determination  of  growth  rate 
quency  reduces  to  the  evaluation  of  volume  integral  over  the 

chamber  and  surface  Integral  over  the  burning  surface.  To  further 
simplify  the  Integration,  we  recall  that 


g(1)  = |jk(u  • VP(1))  - V • (u  * 7q(1)  + q(1)  • Vu)  + V • f(1)] 

(3.25) 

h(1)  = |jkAp(1)  + n • y[o  • 7q(1)  + q(1)  • 7U  - (^)  f (1)]  | (3.31) 


For  the  first  Integral  of  Eq.  3.82,  substituting  Eq.  3.25,  we  have 


g(10)  p(10)*dv 


- I jk(u  • ?P(10))  d/  - 7 • (u  • 7q(10)  + q(10)  • TIj)  p(10)*dv 


+ [ (A)  7 • f (10)  p(10)‘dv 


U * 7p^10^*dv 


(±pao)* 


u • 


dv 


(3.83) 
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Using  Green  s theorem,  the  second  term  can  be  expanded  into  the  form 


/ • p<»V  - { £ 2 • vp<10),2(p(l»)*)dv 


/ <10>Xa  • vp<10))  - (a  • V10))P<10'*j  • n 


ds 


(3.84) 


Recalling  that 


vV10)  = -R2  p(10) 


(3.56) 


110)  . 

n * P - 0 at  the  boundary 


(3.57) 


we  have 


I,  - | (p<10>)*U  • 7p<10>dv  - J-  | [(p<10))*v(u  • Vp<10>)]  • i 


ds 


i / (p<10))%  • f(10)  dv 


(3.85) 


In  similar  fashion  the  second  integral  of  Eq.  3.82  leads  to 


i2  - 1 h<10)  p<"»*d. 


= f + j n • ?(u  • '?p(10>)(p*10'*) 


ds  - 


- 45  - 


n ds 


(3.86) 


■ (u)  l f(10>  PU0>‘- 

s 


Recall  that  A Is  an  admittance  function  with  real  and  Imaginary  parts 
as  shown  In  Eq.  3.26.  When  applied  to  the  burning  surface,  we  can  write 
A as  A = a£  + jA^. 

Substituting  Eqs.  3.85  and  3.86  into  Eq.  3.82,  again  using  Green's 
theorem,  we  get: 


,(n>  + Adi) 


R + H 

E2 

a 


| (p(10))*  U • Vp(10)  dv  - + jaJ) 


r™*™'  ds 


\W  2k^ 


: (10  ) 


7p(10>*dv 


(3.87) 


The  first  term  can  be  written  as 


[ p(10>*  u 


Vp(10)dv 


. E 


(10)*  „(10) 


dv 


l [pa°)%a°) 


U 


n ds 


(3.88) 


This  represents  the  average  convection  of  energy  into  the  system  by 
the  mean  flow,  where  the  (p^10^*  P^10^)/2  term  is  just  a measure  of 
energy  per  unit  volume.  The  second  term, 
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ds 


(3.89) 


-j(Ab + J<){ 


p(10)*  p(10) 


represents  the  work  input  associated  with  the  oscillatory  combustion, 
and  the  last  term  reflects  work  done  on  the  system  by  the  body  force 

or  other  additional  forces. 

Therefore  the  growth  rate  of  the  acoustic  wave  A can  be  expressed 
as: 


t r „ , ,(10)*  (10)  N 

A ' ^ f U • n - ^ I 1 2 ds  + l Ad  <3'9!» 

a s 

where  is  one  of  N damping  factors  including  dissipation  due  to 
aluminum  oxide  particles,  viscous  losses,  energy  lost  due  to  visco- 
elastic motion  of  propellant  grain. 

The  terra  with  the  integral  over  the  burning  surface  can  be 
referred  to  as  since  it  represents  the  source  of  driving  energy 

for  acoustic  waves.  In  this  term  it  includes  the  driving  effects  of 
combustion  as  represented  by  the  admittance  function,  driving  of  wave 
motion  by  the  mean  flow,  and  dissipation  due  to  flow  of  acoustic 
energy  through  the  nozzle.  It  is  important  that  all  effects  which 
drive  oscillation  be  reduced  to  a surface  integral  over  the  combustion 
chamber.  Thus,  once  the  mode  shapes,  for  the  chamber  have  been 

determined,  a simple  integration  is  all  that  is  required  to  determine 
the  driving  effects.  Note  that  l)  ■ n is  simply  the  flow  velocity 
normal  to  the  chamber  surface;  this  quantity  is  readily  calculated  at 
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Che  burning  surface  In  terms  of  propellant  burning  rate  and  Is  zero  on 
Inhibited  surfaces. 

3.4.1  Discussion  of  Acoustic  Gain 

From  Eq.  3.90,  the  gain  expression  of  acoustic  waves  can  be 
written  as: 


a s 


P<10>*  p<10> 

2 


ds 


(3.91) 


2 

with  E defined  as  in  Eq.  3.80. 

a 

To  write  pressure  wave  in  its  fundamental  mode,  i.e.,  a = (0, 
0,  1),  use  Eqs.  3.76  and  3.77  to  ,et 


P 


(10) 

(0,0,1) 


(3.92) 


R,  L and  z are  all  dimensionless  quantities  as  defined  on 
p.  41  shortly  after  Eq.  3.75.  Also  using  the  boundary  conditions 
defined  in  Eq.  3.77,  rewriting  Eq.  3.80  with  m * n * 0,  l = 1,  we 
have 


'gain  ‘ ^ { [-  - " ‘ *b]  cos2  {f  z)  ds  °'n) 

s 

Since  in  cylindrical  coordinates  ds  = rdrdifi,  carrying  out  the  integral, 
we  p,et 


gain 


vR 

L 


(2  ' “ - *h)  cos2  (l?  ') 


(3.94) 
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Now,  if  we  summarize  the  derivation  of  time  dependent  pressure  wave 
expression  for  the  fundamental  mode,  we  get 


(z,t) 


P + € 
0 


(V) 


cos 


(”  z z) cos  (~r) 


(3.95) 


where  t is  a dimensionless  quantity  and  t'  ia  dimensional.  Examining 
Eq.  3.95  we  can  see  that  the  term  • n - A^  in  Eq.  3. 94  plays  a 
very  important  role  in  determining  the  growth  or  decay  of  an  acoustic 
wave  inside  the  combustion  chamber. 

At  the  present  time,  the  data  for  the  calculation  of  A^  at  the 
burning  surface  are  very  limited.  Figure  3.1  shows  some  experimental 
results  by  Angelus  ^ for  doublebased  propellant  and  the  calculated 
curve  by  Culick30  for  hollow  cylindrically  caat  propellant.  Both 
results  are  for  very  high  pressure  compared  to  present  work.  However, 
for  oscillation  frequency  on  the  order  of  kilohertz,  the  value  of  A^_ 
is  more  than  likely  to  be  a small  or  even  negative  number. 


3.4.2  Discussion  of  Acoustic  Loss 

The  decrease  of  acoustic  energy  in  the  combustion  chamber  of 
a rocket  motor  can  be  ascribed  either  to  losses  In  the  medium  or  to 
losses  occurring  at  the  system  boundaries.  Various  losses  and  their 
magnitudes  were  compiled  by  Buff man,  Dehority,  Slates,  and  Price.31 
Table  3.2  shows  these  results  in  terms  of  present  notation. 

From  this  table  it  is  clearly  seen  that  nozzle  loss  is  the  most 
important  one.  Also,  in  Buffman,  et  al.,  it  is  suggested  experimentally 
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real  part  of  burning  surface  admittance 


Frequency  (Hr) 


Fig.  3.1,  Response  of  a burning  propellant  to  steady  pressure 
oscillations.  (After  F.  E.  C.  Culick.) 


Table  3.2.  Mechanisms  contributing  to  acoustic  losses. 


yAd(f)  (SeC-1) 

Losses  in  the  medium  (or  bulk  losses) : 

Viscous  bulk  losses  | 

Heat  conduction  bulk  losses) 
Molecular  relaxation  bulk  losses 
Inhomogeneous  media  bulk  losses 

<0.0014 

<0.6 

'0 

Losses  at  the  system  boundaries  (or  wall  losses): 

Viscous  wall  losses  \ 

Heat  conduction  wall  losses) 

Acoustic  radiation  losses  to  the  walls 
Absorption  losses  to  the  porous  piste 
Transmission  losses  through  the  nozzle 

'5 

'0 

'0 

50  to  600 

that  the  longitudinal  wave  nozzle  loss  term  can  be  written,  in  present 
notation,  as: 


v 


E(001) 


nozzle 


2 


ds 


(1.58)  J (3.96) 


'»  = (0,  0,  1) 


of  which 


nozzle  throat  area 

channel  area  at  nozzle  entrance 


For  a wave  transverse  to  the  chamber  axis,  Culick77  proved 
that  nozzle  admittance  may  be  taken  to  be  identically  zero. 
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Particle  damping  by  aluminum  oxide  smoke  is  also  an  important 
loss  mechanism  and  should  be  included  in  the  calculation.  Culick"7 
adapted  the  work  of  Epstein  and  Carhart32  to  this  problem  and  found: 


2 / 2 \ 2 

/ p\/Ru)\  uR  9p 

I 6o  \f)  2^  f°r  9ma11  Partlcles‘  zf  <K  4T 

\ of  \ o / os 


9p^  wRg 

[Ad|  = ' 6q  for  medium  particles,  <<  <<:  1 

1 ^particles  s o 


2 vw  *•  2 

. R w U)R 

□ g g 

for  large  particles,  1 << 

° 4Vo . o 


where 


(3.97) 


p c v n’ 

r o m o 

50  ‘i  ~J 

s u)R 

s 


p = average  gas  density 


og  * average  solid  density 


C = mass  fraction  of  solid  particles  (or  liquid  particles) 
m 

in  the  combustion  gases 


v = kinematic  viscosity 
o 


= circular  frequency  of  wave  motion  = ^aQ/L^ 


R = average  radius  of  solid  particles 

s 
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Other  loss  chsrscteristics  of  potentlsl  Importance  were  in- 
vestigated during  the  psst.  Hart,  McClure,  and  Bird33  made  a thor- 
ough study  of  the  damping  due  to  the  viscoelastic  motion  of  the  chamber 
boundary.  The  interaction  of  acoustics  with  the  burning  propellant 

surface  was  studied  by  Hart  and  Cantrell.314  Sound  wave  attenuation 

3 5 

in  a doublebased  solid  propellant  was  investigated  by  Nall.  For 
the  present  work,  only  losses  due  to  nozzle  and  particle  damping  will 
be  considered. 
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IV.  HEAT  LOSS  INSIDE  COMBUSTION  CHAMBER 


In  this  chapter,  both  conducting  and  radiating  heat  loas  in- 
side the  combustion  chamber  of  a rocket  motor  vere  investigated  in  a 
very  primitive  form.  Due  to  the  complexity  of  the  combustion  process 
and  mass  and  heat  transfer  phenomena,  several  Important  assumptions 
were  made  to  reduce  the  task  to  a workable  form.  These  assumptions 
include  the  combustion  aystem  as  well  as  some  gas  properties. 

To  theoretically  predict  the  temperature  profile  away  from  the 
burning  surface  of  solid  propellant  inside  the  rocket  combustion  chamber 
is  almost  impossible  without  some  assumptions  about  the  combustion 
system  and  the  properties  of  the  gas.  In  this  atudy,  the  solid  propel- 
lant is  assumed  to  be  burned  uniformly  and  the  burning  aurface  is  always 
normal  to  the  longitudinal  axis  of  the  combustion  chamber.  At  any  cross 
section,  the  gas  temperature  is  assumed  to  be  constant  from  the  center 
to  the  edge  of  the  gas  column.  Conaidering  the  burning  aurface  covera 
the  entire  cylinder,  this  assumption  is  a good  one.  Limited  by  the 
available  data,  some  properties  of  the  gas  have  to  be  assumed.  These 
properties  include  the  flow  velocity,  viscosity,  specific  heat  ratio, 
heat  transfer  coefficient,  and  emissivity  of  the  gas. 

Consider  a heat  flow  system  with  the  fluid  flow  inside  a heat- 
conducting  cylindrical  conduit.  Assume  the  inside  wall  temperature  ia 
room  temperature  before  the  gas  passes  through,  and  the  fluid  temper- 
ature is  much  higher  than  the  room  temperature.  The  heat  losses  when 
fluid  flows  through  the  conduit  are  due  to  conduction  as  well  as  to 


radiation. 


The  conduction  heat  loss  across  the  solid-fluid  interface  is 


dependent  on  the  area  of  the  interface  and  the  temperature  difference 
between  them.  Hence  we  can  write  the  heat  flow  out  of  the  fluid  due 
to  conduction,  Qc,  as 

Qc  = hAAT  (4.1) 

where  h is  the  heat  transfer  coefficient,  A is  the  interface  area, 
and  AT  is  the  temperature  difference.  Rewriting  Eq.  4.1  and  referring 
to  Fig.  4.1,  we  have 


dQc  * h(«Ddz)  (Tf  - Tw)  (4.2) 

Tj  is  the  fluid  temperature  and  Ty  is  the  inside  wall  temperature  of 
the  conduit. 
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The  gas  generating  rate,  G,  of  burning  propellant  is  defined 


G - pAU 


(4.3) 


where  p is  the  density  of  the  solid  propellant,  is  the  burning 
area,  and  U is  the  burning  rate  of  the  propellant. 

As  the  gas  travels  down  the  conduit,  the  heat  transfer  phenom- 
ena can  be  described  as 


GC  dT.  = nDh  (T.  - T \dz 
p f V f v) 


(4.4) 


where  Cp  is  the  specific  heat  at  a constant  pressure.  Simplifying 
Eq.  4.4,  we  get 


dT,  = | . nh,„  (T,  » T 1 
f D [pU]C  f wy 


(4.5) 


For  the  combustion  chamber,  a good  average  value  for  the  term  h/[pU]C 

I 

was  picked, 3f>  and  Eq.  4.5  rewritten  as 


dT,.  = £ 0.023  N*0,2  (T,  - T ) dz 
t D Re  V f vj 


(4.6) 


where  is  the  Reynolds  number  which  can  be  written  as 


Re  u 


(4.7) 


where  u is  the  viscosity  of  the  gas. 

Besides  the  heat  transferred  by  conduction  to  the  rocket  casing 
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from  burned  gas,  heat  also  radiated  Into  the  solid  steel  chamber. 
This  probably  Is  a more  Important  parameter  to  be  considered.  Un- 
fortunately, due  to  the  complicated  gas  composition.  It  Is  almost 
Impossible  for  us  to  estimate  the  total  emisslvity  of  the  gas.  With 
a different  amount  of  seeded  aluminum,  the  emisslvity  should  be  con- 
siderably different.  For  a simple  heat  radiating  case,  as  suggested 
by  most  transport  phenomena  books,37  we  can  write 

Qr  = oeA(AT)4  (4.8) 

where  Qr  Is  the  heat  flow  out  of  gas  due  to  radiation,  o Is  Stefan- 

Boltzmann  constant,  e Is  emisslvity,  and  A Is  the  Interface  area. 

AT  Is  the  temperature  difference  between  fluid  and  wall.  Since  wall 

temperature  Is  considerably  lower  than  fluid  temperature.  It  Is 

4 4 

reasonable  to  let  (AT)  3 T^. 

Writing  Eq.  4.8  In  differential  form,  and  also  referring  to 
Fig.  4.1,  we  have 


GC  dT 
P 


0fTfirD  Hz 


(4.9) 


Theoretically,  we  can  combine  Eqs.  4.6  and  4.9  to  get  the  total 
temperature  change  as  a function  of  distance;  e.g.. 


dTf  - 


5°-023  NRe'2(T£ 


_ \ . oeirD  „4 

- T»)  + ST  Tf 


dz 


(4.10) 


A simple  integration  will  give  us  a linear  relation  between  tempera- 
ture and  distance.  The  linear  relation  is  more  or  less  unrealistic. 
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This  is  due  to  the  assumption  that  the  wall  temperature  stays  con- 
stant and  becomes  a lossy  section  only  when  gas  passes  through. 

Another  way  to  estimate  the  temperature  change  inside  the 
combustion  chamber  is  to  consider  the  temperature  difference  be- 
tween the  chamber  and  nozzle  exit.  Assuming  M is  the  Mach  number  at 
the  nozzle,  the  ratio  of  chamber  temperature  Tf  to  the  nozzle  gas 
temperature,  T , is  derived  as38 

-4=  l + H^M2  (4.11) 

T 

where  v is  the  specific  heat  ratio.  Knowing  the  chamber  length  and 
with  the  assumption  that  heat  loss  is  uniformly  distributed  over  the 
entire  length  of  the  combustion  chamber,  the  heat  loss  due  to  a 
certain  length  section  can  be  roughly  decided. 
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V.  APPLICATION  OF  THE  THEORY 


In  Chapters  II,  III,  and  IV  the  wave  attenuation  equation, 
pressure  wave  oscillation  equation,  and  heat  loss  equation  were  derived 
in  a general  form.  In  this  chapter,  each  of  the  equations  is  applied 
to  the  rocket  motor  designed  for  this  investigation.  The  test  unit  will 
be  discussed  first,  since  all  the  calculations  and  assumptions  are 
based  on  it.  Second,  the  gas  composition  inside  the  rocket  chamber 
will  be  listed.  With  these  data,  the  attenuation  equation  can  be 
calculated.  Following  is  the  calculation  of  acoustic  wave  gain  and 
loss  term.  The  amplitude  and  frequency  for  a specific  mode  will  be 
presented.  Finally,  the  temperature  profile  away  from  the  burning 
surface  inside  the  combustion  chamber  will  be  calculated. 

5 . 1 Test  Unit 

5.1.1  Experimental  Setup 

Figure  5.1  shows  two  photographs  of  the  test  rocket  motor  setup 
used  in  this  investigation.  The  construction  of  the  motor  is  shown 
in  detail  in  Fig.  5.2.  The  system  connection  schematic  is  shown  by 
the  block  diagram  in  Fig.  5.3. 

The  microwave  portion  of  the  test  unit  consisted  of  a klystron 
with  its  power  supply.  The  klystron  is  connected  to  a horn  antenna 
through  an  attenuator  and  frequency  meter.  A dielectric  lens  is  in- 
serted into  the  horn  to  focus  the  microwave  output.  After  passing 
through  the  rocket  motor,  the  K-band  signal  was  then  picked  up  by  the 
antenna  on  the  opposite  side.  Two  separate  channels  were  fed  into  the 


Fig.  5.1.  Photographs  of  the  test  unit. 
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Fig.  5.2.  Detail  drawing  of  a rocket  motor. 


Rocket 

motor 


Fig.  5.3. 


, 


Block  diagram  of  experimental  setup. 
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recorder  to  cover  a rsnge  of  0 to  38  dB  attenuation.  The  recording 
signal  was  then  registered  on  a four-chsnnel  Sanborn  strip  recorder 
together  with  one  channel  of  mean  pressure  and  position  of  the  antenna. 

The  mean  pressure  that  wss  fed  into  the  Ssnborn  recorder  was 
detected  by  a Statun  pressure  gauge  which  was  mounted  on  top  of  the 
rocket  motor.  A Kistler  pressure  gauge  was  mounted  along  with  the 
Statun  pressure  gauge  and  fed  into  a Honeywell  recorder  as  a back-up 
recorder  in  case  of  failure  of  the  Statun  gauge.  These  arrangements 
are  shown  in  Fig.  5.1. 

5.1.2  Dielectric  Lens  Design 

Two  polyethelene  dielectric  lenses  were  designed  to  fit  into 
the  horn  antenna  thst  was  originally  designed  by  Hwang.39  The  lens 
gives  a beam  width  of  4.32  mm  in  rsdlus  at  the  focsl  point  three  inches 
away  from  the  horn  surface.  The  E-field  rsdiation  pattern  is  shown  in 
Fig.  5.4.  The  theory  psrt  of  the  design  was  given  by  Hwang.90  Physical 
dimensions  of  this  lens  are  shown  in  Fig.  5.5. 

5.1.3  Sample  Preparation 

The  propellant  used  in  this  investigation  was  prepared  by  the 
personnel  of  the  Department  of  Chemical  Engineering,  University  of  Utah, 
Salt  Lake  City,  Utah.  The  propellant  was  mixed  and  cast  under  the  fol- 
lowing specifications: 

Mixing  temperature  65°  C 

Mixing  pressure  0.3"  Hg 
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radiation  pattern  of  horn  antenna  with  dielectric  lens 


Mixing  time 


55  minutes 


Casting  pressure 

0. 3”  Hg 

Curing  temperature 

88°  C 

Curing  time 

7 days 

The  chemical  composition  of  different  propellants  uaed  in  the  experi- 
ment is  listed  in  Table  5.1. 

Table  5.1.  Propellant  composition  (weight  percentage). 


Batch 

Coarse1 2 3 

NH.C10. 
4 4 

15  micron 

nh4cio4 

KP2 

PBAA3 

EP0N2824 

Al5 

1 

38.5 

38.5 

5 

15.3 

2.7 

0 

2 

36.0 

36.0 

5 

15.3 

2.7 

5 

3 

26.8 

40.2 

5 

15.3 

2.7 

10 

4 

18.6 

43.4 

5 

15.3 

2.7 

15 

1 — Coarse  ammonium  perchlorate  particle  size  -200  y. 

2 — KP:  potassium  perchlorate  particle  size  -100  y. 

3 — PBAA:  polybutadiene-acrylic  acid  copolymer. 

4 — F.PON282:  epoxy  resin. 

5 — Al:  aluminum  particle  size  -5  y. 

The  propellant  was  cast  in  the  shape  of  a 2"  diameter  cylinder. 
Since  the  shrinkage  from  curing  temperature  to  room  temperature  was 
0.5  percent,41  it  fit  in  the  2"  diameter  rocket  motor.42  Moreover, 
because  a uniform  burning  surface  is  desired,  a thin  layer  of  ailicone 
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grease  (Dow  Corning  Compound  11)  was  coated  on  the  back  and  side  of 
the  grain  to  prevent  burning.  This  also  sealed  the  space  between  the 
propellant  and  the  rocket  motor. 

5 . 2 Plasma  Composition 

The  calculation  of  plasma  composition  is  based  on  a computer 
program  furnished  to  the  University  of  Utah  by  the  United  States  Air 
Force  Rocket  Propulsion  Laboratory  at  Edwards  Air  Force  Base  in 
California.  For  a given  chemical  composition  of  propellant  and  combus- 
tion pressure,  this  program  will  print  out  all  the  equilibrium  constants 
and  chemical  components  of  the  resulting  gases.  For  the  propellant 
composition  listed  in  Table  5.1,  the  gas  composition  is  listed  in 
Appendix  I for  various  chamber  pressures  of  interest.  The  equilibrium 
chamber  temperature  given  by  this  program  is  listed  in  Table  5.2. 

5. 3 Computer  Program  for  Calculating  Wave  Attenuation 

The  solution  of  the  plasma  attenuation  equation  was  carried  out 
on  the  Univac  1108  computer  at  the  University  of  Utah  Computer  Center. 
The  program  listing  is  given  in  Appendix  II.  The  program  listing  is 
divided  into  three  major  parts: 

1.  Initialization  of  constants  representing  a specific 
condition  of  gas  and  electron  temperature. 

2.  Calculation  of  the  mean  velocity  terms  R and  D,  electron 
density,  and  attenuation. 

3.  Printing  the  calculated  attenuation  and  other  pertinent 
parameters  of  the  plasma. 
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Table  5.2.  Combustion  chamber  equilibrium  temperature  predicted. 


Chamber 

Pressure 

psi 

Aluminum  Contained 
5Z  10%  15% 

5 

2569 

2692 

2764 

10 

2639 

2744 

2822 

15 

2663 

2774 

2856 

20 

2680 

2794 

2879 

40 

2717 

2841 

2934 

60 

2737 

2867 

2965 

80 

2750 

2885 

2986 

100 

2760 

2898 

3002 

120 

2767 

2908 

3014 

5.3.1  Initialization  of  Constants 

The  Initialization  section  defines  physical  constants  such  as 
frequency,  velocity  of  light,  etc.,  and  then  reada  from  the  data  deck 
the  mole  fraction  concentrations  of  the  important  constituenta  of  the 
combustion.  Values  used  for  the  representative  mole  fractions  were 
taken  from  a free  energy  analysis  of  the  combustion.  For  the  propel- 
lant compositions  and  chamber  pressure  described  earlier  in  this  re- 
port, more  than  98  percent  of  the  equilibrium  mass  is  accounted  for  by 
six  gaseous  species: 

1.  Diatomic  hydrogen 

2.  Water 
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3.  Diatomic  nitrogen 

4.  Carbon  monoxide 

5.  Carbon  dioxide 

6.  Hydrogen  chloride 

The  mole  fractions  for  these  constituents  were  taken  from  the  program 
printout  and  constitute  the  values  in  the  summation  of  Eq.  2.29, 
which  expression  is  used  in  the  calculation  of  the  mean  velocity  terms 
B and  D. 

After  the  data  are  read  in,  the  temperature  is  initialized  at 
2000  degrees  absolute.  This  temperature  does  not  represent  the  equili- 
brium temperature  of  the  gas  predicted  by  the  free  energy  analysis,  but 
it  merely  serves  as  a starting  reference  temperature  for  attenuation 
calculations.  During  steady-state  conditions  of  combustion,  it  can  be 
assumed  that  the  electrons  and  gas  are  in  thermal  equilibrium,  yet  it 
is  conceivable  that  in  a transient  situation,  the  temperature  could 
change  from  the  equilibrium  value  and  alter  the  attenuation.  For  this 
reason,  a range  of  temperatures  in  a neighborhood  about  the  equilibrium 
gas  temperature  is  used;  hence,  for  a given  gas  state,  a number  of  at- 
tenuations will  be  calculated  for  a group  of  surrounding  temperatures. 
Two  other  chemical  concentrations  are  read  from  the  data  deck.  These 
are  the  equilibrium  amounts  of  potassium  and  chlorine,  the  assumed 
source  and  sink  of  electrons,  respectively. 

5.3.2  Calculation  of  Attenuation  Quantities 

After  the  physical  constants  of  the  plasma  have  been  defined, 
it  is  necessary  to  calculate  three  quantities: 
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1.  The  mean  velocity  terms  B and  D. 

2.  The  number  denaity  of  electrona, 

3.  The  theoretical  attenuation  itself. 

First  is  the  evaluation  of  the  integrala  of  B and  D from  Eqs.  2.41 
and  2.42.  The  integrands  are  continuous,  positive  functions,  and  the 
integrals  are  easily  evaluated  by  a simple  trapezoidal  method. 

The  number  density  of  free  electrons  is  next  determined  from 
Eq.  2.53  using  a subroutine  to  solve  for  the  roots  of  the  cubic  equa- 
tion. It  must  be  ststed  that  for  some  outlying  temperatures,  the 
coefficients  of  the  cubic  equation  were  such  that  round-off  error  in 
the  cubic  calculations  gsve  no  positive  real  roots.  This  fact  is 
equivalent  to  the  obviously  nonphysical  situation  of  assuming  a 
temperature  far  away  from  the  equilibrium  value.  The  relatively  large 
amount  of  chlorine  compared  to  potassium  rendered  the  cubic  equation 
root  locus  sensitive  to  lower  than  equilibrium  temperatures.  The 
occurrence  of  no  positive  solutions  was  not  felt  important  because 
the  temperatures  giving  no  roots  were  alwaya  more  than  one  hundred 
degrees  from  the  indicated  equilibrium  value. 

Early  in  the  experiments  of  this  investigation,  it  was  seen 
that  the  assumption  of  no  electron  sinks  gave  predicted  attenuations 
that  were  much  larger  than  those  observed  in  the  test  motor;  there- 
fore, the  use  of  Saha's  equation  will  be  assumed  incorrect,  and  the 
computer  program  uses  only  the  cubic  expression  to  determine  the 
electron  density. 
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5.3.3  Theoretical  Reaulta 


From  the  calculated  attenuations  of  the  computer  program 
and  the  equilibrium  temperatures  predicted  by  the  Air  Force  program, 
three  plasma  attenuation  characteristics  may  be  derived.  One  set  of 
curves  represents  nonequilibrium  states  of  the  plasma  and  is  a graph 
of  attenuation  versus  temperature  over  a temperature  range  about 
equilibrium.  The  other  two  graphs  are  for  an  assumed  equilibrium 
state,  predicted  by  the  Air  Force  program;  one  shows  attenuation 
versus  temperature  for  the  various  compositions  of  propellant,  and 
the  other  shows  attenuation  versus  pressure  for  different  amounts  of 
aluminum  contained. 

Figures  5.6  through  5.9  are  plots  of  the  nonequilibrium  states 
of  the  plasma  taken  directly  from  the  values  printed  by  the  attenua- 
tion program  for  the  different  concentration  of  ai'iminum.  From  these 
graphs,  it  is  seen  that  the  predicted  attenuation  is  a strong  func- 
tion of  both  the  aluminum  content  and  pressure.  Also  these  graphs 
clearly  show  how  much  attenuation  is  to  be  expected  for  a given 
K-band  signal  incident  upon  a 3.975-inch  cross  section  of  homogeneous 
plasma  of  assumed  composition  at  a certain  pressure. 

Figure  5.10  shows  the  theoretical  attenuation  versus  temperature 
for  the  equilibrium  condition.  From  the  lines  of  constant  pressure 
that  are  shown,  it  is  again  seen  that  the  attenuation  is  a strong 
function  of  pressure  and  aluminum  contained. 

Also  in  Fig.  5.10,  the  pressure  range  lowered  to  5 psi  to  show 
that  wave  attenuation  peaks  up  in  the  vicinity  of  one  atmosphere 
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Fig.  5.6.  Theoretical  attenuation  versus  temperature  and  pressure 
for  a seeding  of  5 percent  potassium  perchlorate  and  0.0 
percent  aluminum. 
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Attenuation , dB  per  3.975-inch  path 


Fig.  5.7.  Theoretical  attenuation  versus  temperature  and  pressure 
for  a seeding  of  5.0  percent  potassium  perchlorate  and 
5.0  percent  aluminum. 
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Fig.  5.9.  Theoretical  attenuation  versus  temperature  and  pressure 
for  a seeding  of  5 percent  potassium  perchlorate  and  15 
percent  aluminum. 
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Attenuation,  d3,  per  3.975-inch  path 
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Fig,  5.10.  Theoretical  attenuation  versus  temperature  for 
equilibrium  condition. 
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pressure  range.  This  is  true  for  all  three  kinds  of  propellants  con- 
sidered. Figure  5.11  shows  this  result  more  clearly  where  wave  attenua- 
tion is  plotted  versus  pressure  at  equilibrium  condition.  Also  from 
the  computer  printout,  we  found  the  electron  drift  velocitv  peaks  up 
at  one  atmosphere  vicinity  which  is  the  direct  result  of  maximum  wave 
attenuat ion. 

The  experimental  results  will  be  presented  in  the  next  chapter. 

5.4  Calculation  of  Acoustic  Cain  and  Loss 

In  this  section  both  acoustic  gain  and  loss  will  be  calculated 
according  to  the  equations  derived  in  Chapter  III.  For  the  two  firings 
picked  for  this  calculation,  both  show  that  the  loss  is  greater  than 
the  gain  when  only  the  fundamental  mode  is  considered  for  pressure  wave 
and  only  nozzle  and  p.  . -le  damping  loss  are  considered  for  loss 
mechanism.  These  calculations  practically  predicted  the  burning  is 
stable . 


5.4.1  Acoustic  Gain 

The  acoustic  gain  expression  for  fundamental  mode  pressure  waves 
was  given  by  Eq.  3.94.  Rewriting  this  equation,  we  have 


A . 
gain 


vR  / , - ,r\  2 / nR  \ 

r v] ' n " \) cos  It  z) 


(3.94) 


Recall  that  R/L  is  a dimensionless  quantity. 

For  the  firing  dated  April  24,  1973,  a 10  percent  aluminum 
propellant,  the  acoustic  gain  is  calculated  and  plotted  as  a function 
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of  a5.  This  is  shown  in  Fig.  5,12.  For  this  firing  the  total  burning 
time  was  3.8  seconds. 

Figure  5.13  shows  the  results  of  the  firing  dated  March  15,  1973, 
a 5 percent  aluminum  propellant.  The  total  burning  time  was  3.3  seconds. 

5.4.2  Acoustic  Loss 

In  this  section,  only  the  nozzle  loss  will  be  considered  since 
the  particle  damping  loss  is  generally  one  order  of  magnitude  less  than 
nozzle  loss.  From  Eq.  3.96,  the  nozzle  loss  is  a function  of  the  square 
of  nozzle  diameter  for  a given  combustion  chamber.  Calculations  are 
made  for  April  24,  1973  and  March  15,  1973  firings  and  the  results  are 
0.1320  and  0.0903,  respectively.  These  results  are  plotted  on  Figs. 

5.12  and  5.13. 

Referring  to  Fig.  5.12,  a stable  system  obtained  experimentally, 
if  A^  < 0 for  a 10  percent  aluminized  propellant,  the  acoustic  loss 
will  exceed  the  acoustic  gain.  Figure  5.14  shows  the  firing  data  of 
June  20,  1972  where  the  high  chamber  pressure  caused  the  failure  of 
the  quartz  window.  This  window  failure  may  have  resulted  from  the 
high  mean  chamber  pressure  and/or  the  acoustic  oscillation  may  have 
caused  the  build  up  of  chamber  pressure.  Figure  5.15  shows  the  acoustic 
gain  and  loss  calculated  according  to  Eqs.  3.94  and  3.96.  The  acoustic 
loss  shown  is  much  less  than  acoustic  gain  even  when  A^  - -0.5. 

5.5  Theoretical  Calculation  of  Temperature 
Change  Inside  the  Combustion  Chamber 

Although  the  heat  loss  equation  has  been  derived  in  Chapter  IV, 


79  - 


5 24  L Nozzle  size  = 0.5781"  diame 


Acoustic  gain  for  April  24,  1973  firing. 
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Fig.  5.13.  Acoustic  gain  for  March  15,  1973  firing. 
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Fig.  5.14.  Pressure  wave  of  firing  dated  June  20,  1972. 


Acoustic  gain  for  June  20,  1972  firing. 


It  is  still  very  difficult  to  carry  out  the  calculation.  The 
main  problem  is  that  there  is  not  enough  information  about  the  gaa 
property  inside  the  combuation  chamber.  These  properties  include  gas 
viscosity,  flow  apeed,  emisaivity,  and  the  thermal  conductivity  of  the 
motor.  Using  some  of  the  "guesaed"  numbers,  Eqa,  4.10  and  4.11  were 
calculated  for  both  ten  and  five  percent  aluminum  propellant. 

For  all  calculations  emissivity  is  chosen  to  be  0.5.  Viscosity 
of  the  gas  is  0.17  lbffl/ft  hr.  For  specific  heat  the  value  given  by  the 
computer  program  from  Edwards  Air  Force  Base  was  adopted.  These  values 
are  listed  in  Appendix  I. 

For  each  kind  of  propellant,  only  one  value  of  gas  density  was 
used  since  the  gas  composition  as  a whole  does  not  vary  too  much  in 
the  pressure  range  of  interest. 

The  results  of  the  calculation  are  listed  in  Table  5.3. 

Table  5.3.  Calculated  temperature  change  inside  the 
combustion  chamber  in  one  inch  distance. 


Type  of 
Propellant 

From  Eq.  4.10 

From  Eq.  4.11 

5%  A1 
10%  A1 

76.5°  C 
77.2°  C 

45.1°  C 
45.1°  C 
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VI.  EXPERIMENTAL  RESULTS 


In  this  chapter,  tie  experimental  results  are  presented  to 
support  the  theoretical  predictions  in  Chapters  II,  III,  and  IV.  First, 
it  will  be  shown  that  the  rocket  chamber  temperature  can  be  measured  by 
microwave  technique  accurately,  even  when  the  solid  propellant  was 
seeded  with  as  much  as  ten  percent  aluminum  by  weight.  Then,  the 
combustion  chamber  temperature  profile  is  plotted  as  a function  of 
distance  relative  to  an  initial  point  away  from  the  burning  surface. 

The  oscillation  of  acoustic  waves  was  observed  by  pressure  wave  detec- 
tion and  microwave  attenuation  measurement,  and  the  oscillation  frequency 
is  compared  with  theoretical  calculation.  Finally,  the  growth  of 
acoustic  waves  is  presented  as  a function  of  burning  surface  admittance. 

6 . 1 Chamber  Temperature  Measurement 

As  mentioned  previously,  the  temperature  measurements  are  based 
on  the  assumption  that  both  electron  sources  and  sinks  are  present  In 
the  plasma.  Since  the  experimental  setup  was  discussed  in  Chapter  V, 
only  the  results  will  be  presented  here. 

Figure  6.1  shows  the  firing  data  of  a 5 percent  aluminum  propel- 
lant with  low  chamber  pressure.  The  length  of  the  propellant  was  1.6 
inches  long,  2 inches  in  diameter,  and  the  nozzle  was  0.484  inch  in 
diameter.  Sample  points  were  taken  and  plotted  as  shown  in  Fig.  6.2. 

The  dashed  line  shown  is  the  equilibrium  temperature  for  various 
pressure  ranges.  From  this  figure,  the  equilibrium  temperature  detected 
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Fig.  6.1.  Experimental  data  for  March  7,  1973  firing. 


Attenuation  per  3.975  inches  path  (dB) 


Fig.  6.2.  Plotted  data  corresponding  to  Fig.  6.1. 
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experimentally  for  5 percent  aluminum  at  5 psi  is  2637°  K.  Figure  6.3 
shows  the  firing  data  of  a ten  percent  propellant  with  the  chamber 
pressure  of  40  pai  in  its  equilibrium  condition.  The  propellant  is  2 
inches  long  and  2 inches  in  diameter.  This  will  be  the  case  unless 
otherwise  noted.  Figure  6.4  shows  the  sample  points  plotted  on  the 
theoretically  calculated  curves.  From  this  figure,  the  equilibrium 
temperature  of  the  combustion  chamber  is  determined  to  be  2770°  K. 

Both  results  are  compared  with  the  values  listed  in  Table  5.2.  This 
is  shown  in  Table  6.1. 


Table  6.1.  Equilibrium  temperature  comparison  for  microwave 
measured  value  and  Air  Force  theoretical  specific 
impulse  (ISP)  computer  program  predicted  value. 


Aluminum 

Contained 

Calculated  Value  Using 
Microwave  Measurement 
Results 

Predicted  Value  by  Air 
Force  ISP  Computer 
Program 

Difference 

5% 

2637°  K 

2569°  K 

2.65% 

10% 

2770°  K 

2841°  K 

2.50% 

Several  attempts  were  made  to  detect  the  chamber  temperature 
propellant,  and  the  results  are  very  poor.  The  main  problem  involved 
is  the  accumulation  of  residue  around  the  nozzle  entrance.  This  ac- 
cumulation shrinks  the  nozzle  passage  considerably,  and  causes  the 
pressure  buildup  ins.de  the  combustion  chamber,  which  results  in  the 
opening  of  the  relief  valve  or  the  breaking  of  the  quartz  window  or 
hoth.  Figure  6.5  shows  the  end  piece  of  the  rocket  motor  after  firing. 
The  propellant  used  was  a 10  percent  one  and  the  nozzle  size  was  5/8 
inch  in  diameter  hefore  firing.  In  order  to  compensate  for  the  clog- 
ging effect,  larger  nozzles  were  used.  But  for  low  pressure  measure- 
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Experimental  result  for  April  24,  1973  firing. 
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Fig.  6.5.  End  piece  of  rocket  motor  after  firing. 

ments,  the  microwave  attenuation  went  beyond  the  capacity  of  present 
equipment.  When  the  pressure  is  extremely  low,  it  is  hard  to  decide 
■list  exactly  what  the  pressure  was.  This  situation  can  be  seen  more 
clearly  from  Fig.  5.11. 

More  experimental  determination  of  the  equilibrium  temperature 
can  also  be  obtained  from  the  data  shown  in  the  next  section. 

6. 2 Temperature  Profile  Measurements 

For  this  part  of  the  experiments,  the  antennas  were  scanned 
simultaneously  across  the  window  during  the  firing  process.  The  micro- 
wave  attenuation  measurements  made  cover  a horizontal  distance  range 
inside  the  combustion  chamber  away  from  the  burning  surface  of  the 
propellant.  During  the  data  reduction  process,  the  position  of  the 
antenna  and  the  location  of  the  burning  surface  are  carefully  calculated 
from  the  burning  rate  of  the  propellant  for  each  datum  point,  so  the 
measured  value  is  the  true  value.  The  first  datum  point  was  chosen  on 
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the  basis  that  the  propellant  was  Ignited  and  the  equilibrium  condi- 
tion is,  or  is  nearly,  reached.  Referring  back  to  Fig,  6,4,  the  num- 
bered data  points  represent  the  time  sequence  of  each  point,  (Later 
this  will  be  indicated  by  arrow  marks  with  only  point  1 identified.) 
Including  the  effects  of  antenna  position  and  the  burning  of  propel- 
lant, the  temperature  profile  was  plotted  as  shown  in  Fig.  6.6. 

Point  1 was  chosen  as  a starting  point  and  all  the  distances  between 
points  from  there  on  were  calculated  relative  to  this  starting  point. 

Theoretically  predicted  results  from  Eqs.  4.10  and  4.11  were 
also  shown  in  Fig.  6.6  for  comparison  purpose.  For  10  percent  propel- 
lant, the  results  are  extremely  close.  However,  it  should  be  noted 
thac  in  the  results  from  Eqs.  4.10  and  4.11,  only  the  slope  of  the 
line  is  important. 

Figure  6.7  shows  other  10  percent  propellant  firing  data.  In 
similar  fashion,  sampling  points  were  plotted  in  Fig.  6.8  and  the 
temperature  profile  is  shown  in  Fig.  6.9.  For  5 percent  propellant, 
the  firing  data  are  shown  in  Figs.  6.10  and  6.11,  and  the  profile  is 
shown  in  Fig.  6.12.  More  data  for  5 percent  aluminum  propellant  are 
shown  in  Figs.  6.13  and  6.14. 

6.3  Oscillation  of  Acoustic  Wave  Measurements 

In  this  section,  the  oscillation  of  acoustic  waves  inside  the 
cylindrical  end  burner  will  be  shown,  and  the  oscillation  frequency 
will  also  be  calculated  accordingly.  The  theoretical  prediction  about 
oscillation  frequency  is  calculated  from  Eq.  3.79  and  listed  in  Table 
6.2  for  several  modes  and  chamber  lengths  for  comparison. 
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Experimental  data  for  March  26,  1^73  firing. 


per  3.975  inches  path  (dB) 


Fig.  6.8.  Plotted  data  corresponding  to  Fig.  6.7. 
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Fig.  6.10.  Experimental  results  for  March  15,  1973  firing 


Attenuation  dB,  per  3.975  inches  path 


2900 


Relative  distance 


Fig.  6.13-  Experimental  data  for  March  5,  1973  firing. 


Firing  date  - March  5,  1973 


Plotted  data  corresponding 


Table  6.2.  Theoretically  predicted  oscillation  fre- 
quency. All  frequencies  are  in  kHz. 

Mode 

'1,  0,  1)  (2,  0,  1)  (0,  1,  1) 

4.2686  4,4877  4.6775 

3.6240  3.8796  4.0977 

3.1651  3.4548  3.6982 

For  this  part  of  the  experiments,  the  combustion  chamber  was 
increased  to  5-3/16  Inches.  For  the  firing  dated  November  27,  1973, 
two  pressure  transducers  were  used  as  before.  One  of  them  was  fed  into 
the  four-channel  Sanborn  recorder,  and  the  other  was  connected  to  a 
Hewlett-Packard  model  3960  Instrumentation  recorder.  The  data  registered 
on  the  Sanborn  is  shown  in  Fig.  6.15.  Notice  the  pressure  wave  satu- 
rated the  recorder.  Figure  6.16  shows  the  pressure  wave  registered  on 
the  tape  recorder  and  played  back  at  the  recording  speed  to  a storage 
scope.  Figure  6.17  through  Fig.  6.21  are  some  segments  of  the  pressure 
wave  shown  in  Fig.  6.16  played  back  at  a speed  sixteen  times  slower  than 
the  recording  speed.  There  is  time  separation  between  pictures.  Figure 
6.19  is  the  more  often  shown  waveform  and  the  oscillation  frequency 
calculated  from  this  picture  is  4.39  kHz. 

For  the  firing  dated  December  8,  1973,  three  pressure  transducers 
were  used.  The  additional  one  was  mounted  at  the  end  of  the  rocket 
motor  opposite  the  nozzle.  The  signal  detected  by  this  transducer  was 
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Fig.  6.15.  Experimental  result  for  November  27,  1973  firing. 


Pressure  (psl) 


5 percent  aluminum 


Nozzle  size: 

0 . 650"  diameter 


Fig.  6.16.  Pressure  wave  of  firing  dated  November  27,  1973. 


Tape  speed: 

15/16  ips 

Horizontal  scale: 
Time  — 5 msec/cm 


Fig.  6.17.  Detail  of  Fig.  ft.) ft  (one  of  five). 


Tape  speed: 
15/16  ips 


Fig.  6.18.  Detail  of  Fig.  6.16  (two  of  five). 
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quency: 4.39  kHz 
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Fig.  6.19.  Detail  of  Fig.  6.16  (three  of  five). 


then  fed  Into  the  tape  recorder  channel  No.  2.  This  Is  shown  in  the 
lower  part  of  Fig.  6.23.  It  is  apparent  that  the  oscillation  was  not 
seen  at  the  end.  Figure  6.22  shows  the  data  obtained  by  the  Sanborn 
recorder.  On  the  upper  half  of  Fig.  6.23  is  the  signal  taken  from  the 
other  transducer  which  was  mounted  on  the  top  of  the  rocket  tiotor  about 
the  center  of  the  window.  The  signal  on  this  picture  was  played  back 
from  the  tape  channel  No.  1 at  a speed  two  and  a half  times  slower  than 
the  original  recording  speed.  Figures  6.24  through  6.32  are  segments 
of  the  signal  registered  on  the  upper  half  of  Fig.  6.23  played  back  at 
sixteen  times  slower  speed.  The  oscillation  frequency  calculated  from 
Fig.  6.28  is  4.272  kHz,  and  the  frequency  calculated  from  Fig.  6.32 
decreased  to  3.66  kHz.  This  could  be  for  one  or  both  of  the  following 
two  reasons: 

1.  The  increase  of  chamber  length. 

2.  Different  mode  dominates. 

For  the  following  two  firings,  both  5 percent  propellant,  two 
pressure  transducers,  one  at  the  center  and  one  at  the  end,  were  used. 
Together  with  the  wave  attenuation  measurements,  all  signals  were  fed 
into  the  HP  tape  recorder.  Figure  6-33  is  the  attenuation  measurement 
of  the  January  18,  1974,  No.  1 firing.  The  burning  time  of  this  and 
the  next  firing  are  considerably  longer  than  all  the  previous  propel- 
lants for  some  unknown  reason.  The  only  explanation  is  that  this  is 
a different  batch. 

Figures  6.34  and  6.35  are  two  waveforms  of  the  microwave  signal 
registered  at  different  times  of  Fig.  6.33  as  indicated.  The  micro- 
wave  signal  remains  in  its  undistorted  original  waveform.  This  shows 
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Fig.  6.22.  Experimental  result  for  December  8,  1973  firing. 


Pressure  (pal) 


5 percent  aluminum 


Nozzle  else: 
0.6500"  dimeter 

Tape  speed: 

15  lpa 

Time  scale: 

0.2  sec/cn 

Volts: 

0.5  V,c* 


f Fig.  6.23.  Pressure  wave  of  firing  dated  December  G,  1973. 
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Fig.  6.2<  . Detail  of  upper  waveform  of  Fig.  6.23 
(one  of  nine). 
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I*pt  speed: 
15/16  ipa 

lime  scale: 
5 asec/ca 


Tape  speed: 
15/16  lps 

Time  scale: 
5 msec /cm 


Tape  speed: 
15/16  ips 

Time  scale: 
5 msec /cm 


Fig.  6.27.  Detail  of  upper  waveform  of  Fig.  6.23 
(four  of  nine). 


Tape  speed: 

15/16  lps 

Time  scale: 

5 msec/cm 

Oscillation  fre- 
quency: 4.272  kHz 


Fig.  6.23.  Detail  of  upper  waveform  of  Fla.  6 21 
(five  of  nine). 


Ill 


Tape  speed: 
15/16  ips 


Tine  scale: 
5 nsec/cn 


Fig.  6.29.  Detail  of  upper  waveform  of  Fig.  6.23 
(six  of  nine). 


Tape  speed: 
15/16  Ips 

Time  scale: 
5 msec /cm 


Fig.  6.30.  Detail  of  upper  waveform  of  Fig.  6.23 
(seven  of  nine). 
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Tape  speed: 
15/16  ips 


Time  scale: 
5 msec/cm 


Fig.  6.31.  Detail  of  upper  waveform  of  Fie.  6.23 
(eight  of  nine). 


Tape  speed: 
15/16  ips 

Time  scale: 
5 msec/sec 

Oscillation 
frequency: 
3.66  kHz 


Fig.  6.32.  Detail  of  upper  waveform  of  Fig.  6 23 
(nine  of  nine). 


Attenuation,  dB 


5 percent  A1  propellant 


Nozzle  size: 
0.640"  diameter 

Tape  speed: 

15  ips 

Time  scale: 

5 sec /cm 


Fig.  6.33.  Attenuation  measurements  of  firing  dated 
January  18,  1974,  No.  1. 


Tape  speed: 
15/16  ips 

Time  scale: 
5 msec /cm 


Fig.  6.34.  Waveform  of  Fig.  6.33. 
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that  pressure  wave  oscillation  does  not  affect  the  microwave  signal. 
Figure  6.36  is  the  pressure  wave  registered  by  the  center  transducer. 

By  slowing  down  the  playback  tape  speed  sixteen  times,  the  oscillation 
of  acoustic  waves  can  be  shown  clearly.  Figure  6.37  shows  the  early 
part  of  the  pressure  wave  and  the  oscillation  frequency  is  3.90  kHz. 
Figures  6.38  shows  the  latter  part  of  the  firing  and  the  oscillation 
frequency  is  4.11  kHz.  Between  these  two  high  pressure  peaks  is  a 
region  of  lower  pressure.  Figure  6.39  tells  the  detail  of  this  region 
where  the  oscillatory  curve  is  the  sixty  cycle  noise.  If  we  take  this 
period  as  equilibrium  situation,  the  temperature  predicted  by  Fig.  5.7 
is  2500  k.  Compared  with  the  value  listed  in  Table  5.2,  the  error  is 
6.12  percent. 

For  the  first  time  pressure  ware  oscillation  was  also  observed 
by  the  transducer  at  the  end.  Figure  6.40  shows  th°  signal  recorded 
on  the  tape  and  Fig.  6.41  shows  the  oscillation  waveform.  Although 
both  magnitudes  are  considerably  lower  than  the  center  one,  the  oscil- 
lation of  acoustic  waves  is  clearly  displayed.  The  oscillation  frequency 
shown  in  Fig.  6.41  is  4.4  kHz  which  is  higher  than  the  center  one  ob- 
served. 

The  results  of  this  firing  are  quite  "different”  from  all  the 
other  firings.  This  is  also  the  first  and  only  time  that  the  oscillation 
of  acoustic  waves  "disappears"  during  the  burning  process.  This  "transi- 
tion" period  takes  approximately  five  seconds,  which  is  longer  than  the 
initial  or  the  final  oscillation  period. 

Figure  6.42  shows  the  attenuation  of  microwave  signals  for  the 


- 115  - 


Tape  speed: 
15/16  lps 


Time  scale : 
5 msec/cm 


? 

I 


Tape  speed: 

15/16  ips 

Time  scale: 

5 msec /cm 

Oscillation  fre- 
quency: 390  kHz 


Fig.  6.37.  Pressure  waveform  of  Fig.  6.36. 


Tape  speed: 

15/16  ips 

Time  scale: 

5 msec/cm 

Oscillation  fre- 
onoucy:  4.11  kHz 


Fig.  6.38.  Pressure  waveform  of  Fie.  6.36 
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Pressure  (psl) 


Tape  speed: 
15/16  lps 

Time  scale: 
5 msec/cm 


Fig.  6.19.  Pressure  waveform  of  F i < > . 6.36. 


Tape  speed: 

15  ips 

Time:  5 sec/cm 

Vertical  scale: 
0.05  V/cm 


Fis-,.  6.40.  Pressure  wave  observed  bv  the  end  transducer 
for  .January  Itt,  1974  , No.  1 firing. 
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Attenuation,  dB 


; 


Tape  speed: 

15/16  ipa 

Tine  scale: 

5 nsec/cn 

Vertical  scale: 
0.01  V/cm 

Oscillation  fre- 
quency: 4.4  kHz 
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Fig.  6.41.  Oscillation  waveform  of  Fig.  6.40. 


5 percent  A1  propellant 


Nozzle  size: 
0.575”  diameter 

Tape  speed: 

15  ips 

Tine  scal>>: 

5 sec /cm 


Fig.  6.42.  Microwave  attenuation  for  the  January  18, 
1974  firing  No.  2. 
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January  18,  1974  firing  No.  2.  At  the  center  part,  the  attenuation  is 
about  4 dB  down.  Figure  6.43  shows  he  puzzling  pressure  wave  registered 
on  the  tape  by  the  center  transducer.  This  could  be  a result  of  the 
malfunctioning  of  the  transducer.  Figure  6.44  is  the  signal  obtained 
by  the  end  one.  The  oscillation  phenomenon  was  again  observed.  This 
is  displayed  in  Fig.  6.45,  and  the  oscillation  frequency  is  4 kHz. 

This  is  very  close  to  the  result  of  the  No.  1 firing.  Assuming  Fig. 

6.43  is  the  true  pressure  wave,  the  oscillation  frequency  obtained 
from  Fig.  6.46,  the  slowing  down  segmer*  of  Fig.  6.41,  is  3.2  kHz. 

From  the  two  firing  results  mentioned  above,  it  is  obvious  that 
at  different  points  of  the  combustion  chamber,  different  modes  of 
acoustic  wave  oscillation  will  be  observed. 


5 percent  A1  propellant 


Nozzle  size: 
0.57'  ' diameter 

Tape  speed: 

15  ips 

Time  scale: 

5 sec/cm 


Fig.  6.4"'.  Pressure  wave  for  January  18,  1974  firing  No.  2. 
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Pressure  (psi) 


i 


* 

t 


f 


Tape  apeed: 

IS  lps 

Time  acale: 

5 sec /cm 

Vertical  scale: 
0.05  V/cm 


Fig.  6.44.  Pressure  wave  recorded  by  the  end  transducer  for 
the  January  18,  1974,  No.  2,  firing. 
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Tape  speed: 
15/16  ips 

Tine  scale: 

5 msec/cra 

Vertical  scale: 
0.005  V/cm 

Oscillation  fre' 
quency:  4 kHz 


Fig.  6.45.  Pressure  oscillation  seen  by  the  end  pressure 
transducer. 
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VII.  SUMMARY  AND  CONCUJSIONS 


On  previous  pages,  the  theory  of  using  microwave  attenuation 
measurement  techniques  as  a plasma  diagnostic  tool  was  discussed.  ?.nd 
the  experimental  results  proved  Its  accuracy  and  effectiveness  even 
when  the  plasma  la  clouded  with  metallic  particles.  The  results  In 
Table  6.1  show  that  the  experimental  values  are  within  three  percent 
of  the  predicted  theoretical  values.  This  means  that  a temperature 
change  of  several  degrees  out  of  several  thousand  degrees  can  be 
measured. 

The  advantages  of  microwave  attenuation  measurement  techniques 
are  as  follows: 

1.  Relatively  simple  equipment  is  needed  as  compared  to  the 
equipment  for  optical  methods. 

2.  Measurements  are  made  directly  on  any  Internal  part  of  the 
plasma.  This  is  done  by  focusing  the  microwave  radiation 
to  a few  millimeters  In  radius,  as  described  in  Chapter  V. 
The  two-color  method  is  presently  the  most  common  rethod 
used  to  measure  the  rocket  combustion  chamber  temperature 
with  suspended  metallic  particles.  This  technique  is  a 
measurement  made  of  the  flame  radiance  without  penetrating 
the  flame.  The  flame  temperature  is  determined  by  comparing 
the  ratio  of  the  measured  spectral  radiances  at  two  wave- 
lengths with  the  ratio  of  the  known  blackbody  spectral 
radiances  at  the  two  wavelengths.  This  technique  also  as- 
sumes that  the  spectral  absorptance  at  the  two  wavelengths 


U equal.  The  temperature  gradient*  inalde  the  flame  are 
not  detectable  by  this  method  and  the  spectral  absorptance 
at  two  wavelengths  la  not  always  th?  same, 

3.  (easurements  are  made  of  the  Instantaneous  and  continuous 
response  of  the  system. 

4.  Results  are  easy  to  Interpret. 

The  conditions  for  the  microwave  attenuation  measurement  tech- 
nique to  be  applicable  are: 

1.  Plane  wave  assumption  Is  justified. 

2.  A known  number  of  free  electrons  are  present  to  cause 
measurable  attenuation. 

3.  The  composition  of  combustion  chamber  gas  Is  known. 

4.  The  chamber  pressure  Is  measured. 

The  temperature  profile  sway  from  the  burning  surface  of  the 
solid  propellant  was  measured  by  the  combination  of  the  scan  of  the 
microwave  signal  across  the  rocket  chamber,  a scan  distance  of  one  Inch, 
snd  by  the  motion  of  the  burning  surface,  a distance  of  one  Inch,  as 
produced  by  the  burning  Itself.  Figures  6.6,  6.9,  6.12,  and  6.14  show 
some  of  the  experimental  results  and  the  calculated  values  obtained 
from  Eqs.  4.10  and  4.11  for  the  temperature  profiles.  The  slopes  have 
an  angular  difference  of  a maximum  of  25  degrees  for  5 percent  aluminized 
propellant  and  18.5  degrees  for  10  percent  aluminized  propellant. 

The  osclllstlon  of  acoustic  waves  Inside  the  combustion  chamber 
Is  theoretically  explored  and  experimentally  observed.  For  a cylindrical 
end  burner,  the  oscillation  phenomenon  Is  displayed  In  Figs.  6.19  to 
6.46,  and  the  oscillation  frequencies  obtained  experimentally  range  from 
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3.2  to  4,4  kHz,  whereas  the  calculated  values  listed  in  Table  6.2 
range  from  2.98  to  5,13  kHz  for  varloua  oscillation  modes. 

The  analysin  of  acoustic  stability  of  a rocket  motor  requires 
the  understanding  of  the  acoustic  gain  and  loss  mechanisms.  The 
acoustic  gain  and  loss  expressions  are  discussed  In  Chapter  111  and 
applied  to  the  cylindrical  rocket  motor  In  Chapter  V.  Figures  5.12 
and  5.13  show  the  application  of  the  theory  to  the  firing  data  for  two 
firings  where  the  systems  are  stable  and  the  acoustic  loaa  exceeded 
the  gain  (/£  < 0)  in  Fig.  5.12  and  (a£  < -0.25)  In  Fig.  5.13.  When 
the  chamber  pressure  was  extremely  high  and  exceeded  the  quartz  window 
strength,  the  windows  failed  aa  a result  of  one  or  both  of  the  following 
possible  conditions: 

1.  The  chamber  mean  pressure  waa  too  high. 

2.  Acoustic  oscillation  caused  the  pressure  buildup. 

Figure  5.15  shows  a case  of  failure  where  the  acoustic  loss  was  much 
less  than  the  acoustic  gain  (A^  ■ -0.5)  and  failure  could  have  resulted 
from  the  acoustic  oscillation. 

This  investigation  shows  that  microwave  attenuation  measurement 
techniques  can  be  used  to  analyze  the  thermodynamic  properties  of  a 
combustion  plasma  containing  auapended  metallic  particles  in  a solid 
propellant  rocket  motor.  The  acoustic  oscillation  phenomenon  of  an 
end  burner  waa  also  investigated  theoretically  and  experimentally. 
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APPENDIX  I 


COMBUSTION  CHAMBER  GAS  PROPERTIES 


Propellant  A1 
Contained 


0 Percent 


5 Percent 


10  Percent 


15  Percent 


Density 

■ ■■  - ■ 

1.714 

1.735 

1.756 

Heat  Cap.  (cal.) 

43.526 

43.732 

43.770 

Entropy  (cal.) 

280.692 

275.736 

269.927 

Molecular 

Weight 

24.325 

24.836 

25.346 

Total  Moles 
of  Gas 

4.111 

4.026 

3.945 

Moles  of  H2 

0.82705 

1.07610 

1.32189 

Moles  of  HjO 

0.93208 

0.58817 

0.25539 

Moles  of  N2 

0.32037 

0.29922 

0.27824 

Moles  of  C02 

0.18379 

0.09346 

0.03393 

Moles  of  HC1 

0.55497 

0.49985 

0.42980 

Moles  of  CO 

1.03658 

1.12690 

1.18641 

Moles  of  K 

0.00155 

0.00267 

0.00417 

Moles  of  Cl 

0.05649 

0.06403 

0.06381 

Moles  of  A1 

0.00001 

0.00011 

0.00109 

Moles  of  AljOj 

0.09133 

0.18119 

0.26014 

COMBUSTION  CHAMBER  PRESSURE  TO  PSI 


Contained 


Density 

Heat  Cap.  (cal.) 

Entropy  (cal.) 

Molecular 

Weight 

Total  Moles 
of  Gas 

Moles  of  L 


Moles  of  HjO 


Moles  of  N- 


Moles  of  C0r 


Moles  of  HCl 


Moles  of  CO 


Moles  of  K 


Moles  of  Cl 


Moles  of  A1 


Moles  of  A^O^ 


5 Percent 

10  Percent 

15  Percent 

1,714 

1.735 

1.756 

43.589 

43.787 

43.811 

275.042 

270.205 

264.509 

24.426 

24.966 

25.502 

4.094 

4.006 

3.921 

0.82727 

1.08163 

1.33370 

0.93898 

0.59376 

0.25956 

0.32042 

0.29925 

0.27825 

0.18229 

0.09228 

0.03356 

0.56305 

0.50751 

0.43492 

1.03808 

1.12806 

1.18677 

0.00128 

0.00225 

0.00359 

0.04797 

0.05533 

0.05545 

0.00001 

0.00009 

0.00099 

0.09132 

0.18105 

0.25950 
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Propellant  A1 
Contained 


0 Percent 


5 Percent 


10  Percent 


15  Percent 


Density 

1.714 

1.735 

1.756 

Heat  Cap.  (cal.) 

43.623 

43.817 

43.835 

Entropy  (cal.) 

271.747 

266.982 

261.355 

Molecular 

Weight 

24.482 

25.039 

25.591 

Total  Moles 
of  Gas 

1 

4.085 

3.994 

3.908 

Moles  of  Hj 

0.82732 

1.08466 

1.34046 

Moles  of  H20 

0.94284 

0.59696 

0.26180 

Moles  of  N2 

0.32046 

0.29927 

0.27825 

Moles  of  CO^ 

0.18153 

0.09166 

0.03334 

Moles  of  HCl 

0.56755 

0.51183 

0.43783 

Moles  of  CO 

1.03883  1 

1.12867 

1.18697 

Moles  of  K 

1 

0.00113  ! 

0. 00202 

0.00327 

Moles  of  Cl 

0.04326 

0.05044 

0.05077 

Moles  of  A1 

0. 00001 

0.00008 

0.00092 

Moles  of  AljOj 

0.09133 

| 0.18101 

0.25923 

COMBUSTION  CHAMBER  PRESSURE  20  PSI 


Propellant  A1 
Contained 

0 Percent 

5 Percent 

10  Percent 

15  Percent 

Density 

1.694 

1.714 

1.735 

1.756 

Heat  Cap.  (cal. ) 

43.405 

43.646 

43.839 

43.852 

Entropy  (cal.) 

273.398 

269.413 

264.701 

259.124 

Molecular 

Weight 

23.936 

24.520 

25.090 

1 

25.654 

Total  Moles 
of  Gas 

4.178 

4.078 

3.986 

3.898 

Moles  of  Hj 

0.59052 

0.82732 

1.08670 

1.34517 

Moles  of  H^O 

1.26316 

0.94552 

0.59918 

0.26328 

Moles  of 

0.34183 

0.32049 

0.29929 

0.27826 

Moles  of  CO^ 

0.31150 

0.18105 

0.09126 

0.03318 

Moles  of  HC1 

0.62446 

0.57062 

0.51482 

0.43984 

Moles  of  CO 

0.90887 

1.03930 

1.12906 

1.18711 

Moles  of  K 

0.00051 

0.00103 

0.00187 

0.00305 

Moles  of  Cl 

0.03156 

0.04004 

0.04707 

0.04755 

Moles  of  A1 

— 

0.00001 

0.00008 

0.00086 

Moles  of  A^O^ 

0.09133 

0. 18099 

0.25909 
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COMBUSTION  CHAMBER  PRESSURE  40  PSI 


Propellant  A1 
Contained 

0 Percent 

" ■ - — 
5 Percent 

10  Percent 

15  Percent 

Density 

1.694 

1.714 

1.735 

1.756 

Heat  Cap.  (cal.) 

43.454 

43.698 

43.888 

43.894 

Entropy  (cal. ) 

267.650 

263.806 

259.224 

253.770 

Molecular 

Weight 

23.995 

24.607 

25.208 

25.801 

Total  Moles 
of  Gas. 

4.168 

4.064 

3.967 

3.876 

Moles  of 

0.58803 

0.82711 

1.09125 

1.35617 

Moles  of  H2O 

1.26862 

0.95165 

0.60433 

0.26648 

Moles  of  Nj 

0.34191 

0.32056 

0.29934 

0.27828 

Moles  of  CO2 

0.31092 

0. 18006 

0.09039 

0.03281 

Moles  of  HC1 

0.63072 

0.57755 

0.52169 

0.44449 

Moles  of  CO 

0.90942 

1.04026 

1.12989 

1.18743 

Moles  of  K 

0.00039 

0.00082 

0.00153 

0.00255 

Moles  of  Cl 

0.02518 

0.03281 

0.03938 

0.04019 

Moles  of  A1 

— 

0.00001 

0. 00006 

0.00072 

Moles  of  A^O^ 

— 

0.09136 

0.18098 

0.25893 
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COMBUSTION  CHAMBER  PRESSURE  60  PSI 


Propellant  A1 
Contained 

0 Percent 

5 Percent 

10  Percent 

15  Percent 

Density 

1.694 

1.714 

1.735 

1.756 

Heat  Cap.  (cal.) 

43.479 

43.725 

43.915 

43.919 

Entropy  (cal.) 

264.294 

260.534 

256.032 

250.652 

Molecular 

Weight 

24.026 

24.654 

25.273 

25.884 

Total  Moles 
of  Gas 

4.162 

4.056 

3.957 

3.863 

Moles  of  H2 

0.58669 

0.82689 

1.09365 

1.36234 

Moles  of  H20 

1.27148 

0.95499 

0.60718 

0.26810 

Moles  of 

C. 34196 

0.32061 

0.29937 

0.27830 

Moles  of  C02 

0.31069 

0.17959 

0.08995 

0.03253 

Moles  of  HCl 

0.63398 

0.58126 

0.52544 

0.44705 

Moles  of  CO 

0.90965 

1.04071 

1.13030 

1.18761 

Moles  of  K 

0.00034 

0.00072 

0.00135 

0.00229 

Moles  of  Cl 

0.02187 

0.02894 

0.03519 

0.03616 

Moles  of  A1 

— 

— 

0.00005 

0.00064 

Moles  of  A1203 

— 

0.09138 

— 

0.18101 

0.25894 
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COMBUSTION  CHAMBER  PRESSURE  80  PSI 


Propellant  A1 
Contained 

0 Percent 

5 Percent 

10  Percent 

15  Percent 

Density 

1.694 

1.714 

1.735 

1.756 

Heat  Cap.  (cal.) 

43.496 

43.745 

43.933 

43.936 

Entropy  (cal. ) 

261.916 

258.217 

253.771 

248.445 

Molecular 

Weight 

24.046 

24.685 

25.318 

25.941 

Total  Moles 
of  Gas 

4.159 

4.051 

3.950 

3.855 

Moles  of  H2 

0.58577 

0.82669 

1.09521 

1.36657 

Moles  of  HjO 

1.27336 

0.95724 

0.60914 

0.26913 

Moles  of  N2 

0.34199 

0. 32064 

0.29940 

0.27830 

Moles  of  CO ^ 

0.31056 

0.17930 

0.08967 

0.03243 

Moles  of  HC1 

0.63610 

0.58372 

0.52798 

0.44878 

Moles  of  CO 

0.90977 

1.04098 

1.13056 

1.18772 

Moles  of  K 

0.00030 

0.00065 

0.00123 

0.00211 

Moles  of  Cl 

0.01972 

0.02637 

0.03236 

0.03344 

Moles  of  A1 

— 

— 

0.00005 

0.00058 

Moles  of  AljO 

0.09140 

0.18104 

0.25899 
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COMBUSTION  CHAMBER  PRESSURE  100  PSI 


P’opellant  A1 
Contained 

0 Percent 

5 Percent 

10  Percent 

15  Percent 

Density 

1.694 

1.714 

1.735 

1.756 

Heat  Cap.  (cal.) 

43.507 

43.759 

43.947 

43.950 

Entropy  (cal.) 

260.072 

256.421 

252.021 

246.738 

Molecular 

Weight 

24.060 

24.709 

25.351 

25.985 

Total  Moles 
of  Gas 

4.156 

4.047 

3.945 

3.848 

Moles  of  H2 

1.27474 

0.95892 

0161060 

0.26986 

Moles  of 

0.34202 

0.32066 

0.29942 

0.27830 

Moles  of  CO2 

0.31049 

0.17911 

0.08947 

0.03231 

Moles  of  HC1 

0.63764 

0.58554 

0.52987 

0.45007 

Moles  of  CO 

0.90983 

1.04116 

1.13074 

1.18781 

Moles  of  K 

0.00028 

0.00060 

0.00114 

0.00197 

Moles  of  Cl 

0.01815 

0.02448 

0.03026 

0.03140 

Moles  0.'  A1 

— 

— 

0.00004 

0.00054 

Moles  of 

— 

0.09141 

0. 18106 

0.25906 

- 134  - 


I 


COMBUSTION  CHAMBER  PRESSURE  120  PSl 


Propellant  A1 
Contained 

0 Percent 

5 Percent 

10  Percent 

15  Percent 

Density 

1.694 

1.714 

1.735 

1.756 

Heat  Cap.  (cal.) 

43.517 

43.769 

43.958 

43.961 

Entropy  (cal.) 

258.567 

254.956 

250.593 

245.344 

Molecular 

Weight 

24.072 

24.727 

25.377 

26.020 

Total  Moles 
of  Cas 

4.154 

4.044 

3.941 

3.843 

Moles  of  H2 

0.58451 

0.82633 

1.09723 

1.37231 

Moles  of  H20 

1.27581 

0.96025 

0.61176 

0.27042 

Molea  of  N2 

0. 34204 

0.32068 

0.29943 

0.27830 

Mole a of  C02 

0.31044 

0.17896 

0.08931 

0.03221 

Moles  of  HC1 

0.63883 

0.58695 

0.53135 

0.45109 

Moles  of  CO 

0.90986 

1.04128 

1.13088 

1.18788 

Molea  of  K 

0.00026 

0.00056 

0.00108 

0.00187 

Moles  of  Cl 

0.01694 

! 

0.02301 

0.02860 

0.02980 

Moles  of  A1 



— 

0.00004 

0.00050 

Molea  of  Al^ 

— 

0.09142 

0.18109 

0.25193 

- 135  - 


APPENDIX  II 


ATTENUATION  CALCULATION  COMPUTER  PROGRAM 
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APPENDIX  III 


COMPARISON  OF  WAVE  TRANSMISSIONS  BETWEEN  HORN  ANTENNAS 
WITH  DIRECT  WAVEGUIDE  CONNECTIONS 

Measurements  of  the  transmission  characteristics  of  the  horn 

* 

antenna  were  made  and  the  results  are  given  in  the  following  table 
for  the  various  cases: 

1.  Free  space  between  antennas 

2.  Transmission  through  rocket  motor  with  quartz  window 

3.  Rocket  motor  in  place  but  with  quartz  window  blocked 


Without 

Dielectric 

Lens 

With 

Dielectric 

Lens 

Difference 

Free  space  between 
antennas 

- 3.7  dB 

- 0.6  dB 

3.1  dB 

Through  rocket  motor 
with  quartz  window 

-10.6  dB 

- 1.6  dB 

9.0  dB 

Rocket  motor  in  place 
but  with  quartz  windows 
blocked 

-19.0  dB 

-25.5  dB 

-6.5  dB 

The  insertion  loss  was  compared  to  a 
and  the  measurements  of  horn  without 


direct  connection  of  waveguide, 
lenses  and  with  lenses  were  made . 


* 


H.  S.  Hwang,  C.  H.  Durney,  and  R.  W.  Grow,  "Final  Report  on  the  Study 
of  Dielectric  Insert  Horn  Antennas",  unpublished  memorandum. 


It  is  apparent  that  the  Insertion  loss  of  the  lens~fitted  horn 
antennas  was  0.6  dB  without  the  motor  in  place  and  1,6  dB  with  the 
motor  in  place  with  transmission  through  the  quartz  windows  of  the 
motor. 
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